Manuscript prepared for J. Name 

with version 4.2 of the LTgX class copernicus.cls. 

Date: 25 March 2013 



m 
o 



(N ■ 
(N ; 

5?: 
4d ■ 

i 

O 

CO 



Recent variability of the solar spectral irradiance and its impact on 
climate modelling 

I. Ermolli 1 , K. Matthes 2 , T. Dudok de Wit 3 , N. A. Krivova 4 , K. Tourpali 5 , M. Weber 6 , Y. C. Unruh 7 , L. Gray 8 , 
U. Langematz 9 , P. Pilewskie 10 , E. Rozanov 11 12 , W. Schmutz 11 , A. Shapiro 11 , S. K. Solanki 4 13 , and T. N. Woods 10 

'INAF, Osservatorio Astronomico di Roma, Monte Porzio Catone, Italy 
2 GEOMAR I Helmholtz-Zentrum fur Ozeanforschung Kiel, Kiel, Germany 
3 LPC2E, CNRS and University of Orleans, Orleans, France 

4 Max-Planck-Institut fur Sonnensystemforschung, 37191 Katlenburg-Lindau, Germany 
5 Laboratory of Atmospheric Physics, Aristotle University of Thessaloniki, Greece 
6 Institut fur Umweltphysik, Universitat Bremen FBI, Bremen, Germany 
7 Astrophysics Group, Blackett Laboratory, Imperial College London, SW7 2AZ, UK 

8 Centre for Atmospheric Sciences, Dept. of Atmospheric, Oceanic and Planetary Physics, University of Oxford, UK 
9 Institut fur Meteorologie, Freie Universitat Berlin, Berlin, Germany 

10 University of Colorado, Laboratory for Atmospheric and Space Physics, Boulder, CO, USA 

1 1 Physikalisch-Meteorologisches Observatorium, World Radiation Center, Davos Dorf, Switzerland 

12 IAC ETH, Zurich, Switzerland 

13 School of Space Research, Kyung Hee University, Yongin, Gyeonggi 46-701, Republic of Korea 



> 

in 

cn 
o 
m 

10 

X 



Abstract. The lack of long and reliable time series of so- 
lar spectral irradiance (SSI) measurements makes an accu- 30 
rate quantification of solar contributions to recent climate 
change difficult. Whereas earlier SSI observations and mod- 
els provided a qualitatively consistent picture of the SSI vari- 
ability, recent measurements by the SORCE satellite sug- 
gest a significantly stronger variability in the ultraviolet (UV) 35 
spectral range and changes in the visible and near-infrared 
(NIR) bands in anti-phase with the solar cycle. A number 
of recent chemistry-climate model (CCM) simulations have 
shown that this might have significant implications on the 
Earth's atmosphere. Motivated by these results, we summa- 40 
rize here our current knowledge of SSI variability and its im- 
pact on Earth's climate. 

We present a detailed overview of existing SSI measure- 
ments and provide thorough comparison of models available 
to date. SSI changes influence the Earth's atmosphere, both 45 
directly, through changes in shortwave (SW) heating and 
therefore, temperature and ozone distributions in the strato- 
sphere, and indirectly, through dynamical feedbacks. We in- 
vestigate these direct and indirect effects using several state- 
of-the art CCM simulations forced with measured and mod- 50 
eled SSI changes. A unique asset of this study is the use of a 
common comprehensive approach for an issue that is usually 
addressed separately by different communities. 

We show that the SORCE measurements are difficult to 
reconcile with earlier observations and with SSI models. Of 
the five SSI models discussed here, specifically NRLSSI, 
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SATIRE-S, COSI, SRPM, and OAR, only one shows a be- 
haviour of the U V and visible irradiance qualitatively resem- 
bling that of the recent SORCE measurements. However, 
the integral of the SSI computed with this model over the 
entire spectral range does not reproduce the measured cycli- 
cal changes of the total solar irradiance, which is an essen- 
tial requisite for realistic evaluations of solar effects on the 
Earth's climate in CCMs. 

We show that within the range provided by the recent SSI 
observations and semi-empirical models discussed here, the 
NRLSSI model and SORCE observations represent the lower 
and upper limits in the magnitude of the SSI solar cycle vari- 
ation. 

The results of the CCM simulations, forced with the SSI 
solar cycle variations estimated from the NRLSSI model and 
from SORCE measurements, show that the direct solar re- 
sponse in the stratosphere is larger for the SORCE than for 
the NRLSSI data. Correspondingly, larger UV forcing also 
leads to a larger surface response. 

Finally, we discuss the reliability of the available data and 
we propose additional coordinated work, first to build com- 
posite SSI datasets out of scattered observations and to refine 
current SSI models, and second, to run coordinated CCM ex- 
periments. 



1 Introduction 

The question of whether - and to what extent - the Earth's 
climate is influenced by solar variability remains central to 
the understanding of anthropogenic climate change. Ac- 
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SORCE & TIMED, 22 April 2004 to 23 July 2010 




Fig. 1. The solar spectral irradiance as inferred from SORCE and 
TIMED observations only, from 22 April 2004 till 23 July 2010. 
(A) shows the average solar spectral irradiance for that period. 
A black-body model has been used to extend the SSI for wave- 90 
lengths beyond 1580 nm. (B) displays the characteristic altitude of 
absorption in the Earth's atmosphere for each wavelength, defined 
as the altitude at which the optical depth equals one. (C) shows the 
relative variability (peak to peak/average) for solar cycle variations 
inferred from measurements obtained between 22 April 2004 and 95 
23 July 2010. Spectral regions, where the variability is in phase 
with the solar cycle (represented by, e.g. the sunspot number or 
the TSI) are marked in red, while blue denotes ranges where the 
variability measured by SORCE is out-of-phase with the solar cy- 
cle. These phases, as well as the magnitude of the variability in the 
UV, are not all reproduced by models and other observations (see 
Sect. [3] as well as Figs. [2] [7] and [8), and thus should be considered 
with care. (D) shows the absolute variability, which peaks strongly 
in the near-UV. 



cording to the 4th assessment report of the Intergovernmen- 
tal Panel on Climate Change, solar variability represents 
abo ut 8% of recen t total net radiative anthropogenic forc- 



6o ing dSolomon et all 120071) . However, there is a large uncer- 



tainty in this figure because several aspects of solar forcing 
and the different mechanisms by which solar variability in- 
fluences the Earth's environment are sti l l poorly understood 



(see e.g.lPap et al l l2004t ICalisesi et all 120071: lHaigh. 2007 



Gray et all |20 1 oTlLockwoodl |20 1 2i and references therein) 



For these reasons, the quantification of solar contribution to 
climate change remains incomplete. This is further high- 
lighted by some of the most recent investigations of solar 
spectral irradiance (SSI) variations and estimates of their 
influence on the Earth's atmosphere based on chemistry- 
climate model (CCM) simulations. 

Regular space-based measurements of the solar irradiance 
started in 1978. The total solar irradiance (TSI), i.e. the 
spectrally integrated radiative power density of the Sun in- 
cident at the top of Earth's atmosphere, has been moni- 
tored almost continuously and wa s found to vary on different 



198U Frohlich and Leanl 12004 



time sca l es (jWillson et al 
iFrohlichi 120091) . Most noticeable is the rj 0.1 % modulation 
of TSI in phase with the 1 1-yr solar cycle. Changes 2-3 times 
larger than this are observed on time scales shorter than few 
days. Measurements of SSI, however, are not continuous 
over the satellite era and until recently have concentrated on 
the ultraviolet (UV) radiation, because of the larger relative 
variability of SSI below 400 nm (Fig.Q]) and the impact of 
these wavelengths on the terrestrial atmosphere through ra- 
diative heating and ozone photochemistry. 

SSI variations differ from those observed in the TSI. The 
variability of visible and NIR bands barely exceeds 0.5 % 
over a solar cycle; in the near UV and shorter wavelengths 
variability increases with decreasing wavelength, reaching 
several percent at 200-250 nm, and seve ral tens of percent, 



and even more, below about 200 nm (e.g. lFlovd et all 12003 



and references therein). These bands are almost completely 
absorbed in the Earth's middle and upper atmosphere (Fig.[T] 
and are the primary agent affecting heating, photochemistry, 
and therefore, the dynamics of the Earth's atmosphere. Vari- 
ations of the solar UV radiation between 120 and 350 nm 
lead to changes in stratospheric ozone and heating that am- 
plify the effect of the UV radiation in the Earth's atmosphere, 
possibly also throu gh indirect mecha nisms (e.g. the "top- 



down" mechanism, Gray et al. , 2010l) . Hence, although the 



UV radiation shortward of 400 nm represents less than 8 % 
of the TSI, its variability may have a significant impact on 
climate. In contrast, the visible and IR bands, which have 
the largest contribution to the TSI, small variations over the 
solar cycle, and no absorption in the atmosphere but in some 
well-defined IR bands, directly heat the Earth's surface and 
the lower atmosphere. The large amount of solar flux at the 
visible and IR bands implies that small flux differences may 
induce important terrestrial consequences. The impact of the 
variability of these bands on the Earth's climate is expected 
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Fig. 2. Relative contribution of the UV (200-400 nm), visible (400-700 nm), near-IR (700-1000 nm) and IR (1000-2430 nm) ranges to 
the TSI change over the solar cycle as derived from measurements and models described in Sects. [2] and [3] For SORCE/SIM, only the 
period between 2004 and 2008 can be considered. For other data and models, the plotted relativ e differences are between so lar maximum an 
minimum. Within e ach waveleng t h bin, f rom left to right: WR-2002 (light green, only in UVl lWoods and Rottmanl [20021) , UARS/SUSIM 
(green, only in UV; iMorrill eTall 1201 lab . NRLSSI (black ), SATIRE-S (blue), COSI (purple), OAR (light blue), SCIAMACHY (brown), 
SORCE (red), and SORCE re-analysis (orange, only in UV; IWoodsLl2012h . The exact wavelength ranges used for SUSIM and SCIAMACHY 
in the UV are 150-400 nm and 240-400 nm, respectively. The possi ble related correction s are, however, expected to li e within 2-3 %. Note 
that fo r the SCIAMACHY-based model, the original values listed bv lPagaran et al.l d2009f) are shown. As discussed by|Kr ivova and Solankil 
j2013h . these values should most likely be corrected by a factor of roughly 1.2. 



to be small, although it may involve amp lification mech 
anisms (e.g. the "bottom -up" mechanism, Ivan Loon et al 



2007tlMeehl et al.U2009l) . 



Early satellite measurements of the solar UV vari- 
ability have shown a qu alitatively consistent behaviour 
( DeLand and Cebula , 20121). which is fairly well reproduced , 



by SSI models (e. g 



Lean et al.. 1997: Krivova et all 12006; 



Unruh et al. , 2012 ; Lean and DeLand , 2012b . This situation 



changed with the la unch of the Spectral Irradiance Moni- 
tor instrument (SIM, Harder et al.l 12005 al) onboard the So- 
lar Radiation a nd Climate Experiment satellite (SORCE, 
Rottmanl 120051) in 2003, which was shortly after the most 1 
recent maximum of solar activity. The SORCE/SIM data 
showed a four to six times greater decrease of the UV radi- 
ation between 200 a nd 400 nm over the period 2004-2008 
dHarder et al.l I2009I) . part of the declining phase of solar 



cycle 23, compared to earlier measurements and models 
dBall et all l201ll IPagaran et al.l |201 lal lUnnih et all I2012; 
DeLand and Cebulal l2012t " Lean and DeLandl I2012 ). This 
larger decrease measured in the UV (Fig. |2}, which exceeds 
the TSI decrease over the same period by almost a factor of 
two, is compensated by an increase in the visible and NIR 



bands. Variability out-of-phase with solar activity is indeed 
predicted by some SSI models in the NIR, but with a sig- 
nificantly lower magnitude than found by SORCE/SIM. De- 
tails are provided in the following. The inverse variability 
observed by SIM in a wide integrated band in the visible 
was, however, unexpected. It can be interpreted as a result 
of effects induced by the evo lution of sur f ace m agnetism 
in the solar atmosphere (e.g. Harder et al.l 2009 ). How- 
ever, other observations and analyses of existing long-term 
SSI data show results in contrast with those derived from 



SORCE/SIM (Wehrli et al, 2012) 



When used as solar input to CCM simulations, 
SORCE/SIM observations lead to significantly larger short- 
wave (SW) heating rates in the upper stratosphere com- 
pared to results obtai ned by using the comm o nly u tilised 



NRLSSI model data dLean et al.l 119971: iLeanl |2000f). and 



a decrease of stratosphe ric ozone above an altitude of 45 km 



during solar maximum dHaigh et al.l bOlfj . These changes 



in radiative heating and ozone photochemistry in the strato- 
sphere also impact the responses of the "top-down" solar 
UV mechanism in th e Earth 's atmosphere and at the surface 
dKodera and KurodaL 120021) . which may depart from cur- 
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rent understa n ding (ICahalan et al. L 2010 L Ineson et"all 201 1 



20121) . 



Merkel et all 1201 it bberlander et all 120121 ISwartz et al 



Validation of the results of these model simula- 



tions with ozone, zonal wind, or temperature measurements 
is difficult, because the data are sparse and do not cover 210 
enough solar cycles. Although some ozone observations 



seem to agree with mod el calculations (e.g. Hai gh et al 



2010t IMerkel et all l201ll) , it should be noted that the SIM 



measurements employed for the analysis covered less than 
one solar cycle and required extra polation over a full cycle, 
and therefore, added uncertainty (IGarciai I2010I) . Also the 215 
transition altitude from in-phase (lower and middle strato- 
sphere) to out-of-phase (upper stratosphere) ozone signals 
with the solar cycle is not consistent among the different 
models and requires further investigations. 

Unfortunately, observations of the full solar spectrum will 22 » 
likely have a multi-year gap before the next generation SSI 
instrument is launched. Based on data presently available, 
a thorough understanding of the impact of SSI on climate 
requires verification and validation of existing SSI measure- 
ments for internal consistency, calculations of middle atmo- 225 
sphere climate models with different reliable scenarios of SSI 
variations, and comparison of measurements and model re- 
sults with climate records, i.e. a study involving the coordi- 
nated work of various research communities, which is part of 
the COST Action ESI 0050 to which most of the authors of 
this paper belong. 

This paper summarises and compares, for the first time, 230 
a large number of SSI observations and models, and dis- 
cusses the impact of these data on Ea rth's climate. A re view 
in this area was recently published bv lGrav et al.l (120101) . but 
their focus was more on different possible forcing mecha- 
nisms. Although the Sun affects the climate system in nu- 235 
merous ways, we here focus on radiative forcing only, with 
particular attention given to the role of the SSI rather than 
that of the TSI, which is still the sole solar input in many cli- 
mate models. Since space-based observations of the SSI and 
of the terrestrial atmosphere conditions are sparse or absent 240 
before 1980, we restrict our analysis to the data of the last 
three decades, i.e. roughly the last three solar cycles. Fur- 
thermore, we concentrate on the effects of the UV variability 
because of its potentially large impact on the terrestrial at- 
mosphere. 245 

The paper is organised as follows: in Sect. 2, measure- 
ments of SSI variations are described, and their accuracy 
on time scales from days to 11 -yr solar cycle is discussed. 
In Sect. 3 we delineate mechanisms responsible for the SSI 
variations, outline methods of irradiance reconstructions and 250 
briefly describe and compare several of the most broadly 
used models of SSI variations. Section 4 discusses the impact 
of the current lower and upper boundaries of SSI solar cycle 
estimated variations on the atmospheric response in several 



current state-of-the art CCMs. A summary and concluding 
remarks are provided in Sect. 5. 



2 Solar irradiance measurements 

We first discuss the available solar irradiance measurements, 
focussing on the most recent data. These data are input to 
and, at the same time, the main source of constraints for 
both the SSI models and CCM simulations presented below. 
We describe the evolution of measurements carried out from 
space since 1978 and discuss major instrumental and mea- 
surement problems that limit the creation of single composite 
time series from existing records (Sect. |2~TV Long and accu- 
rate SSI time series are critical to obtain reliable estimates 
of the solar variability impact on the terrestrial atmosphere. 
We then present instruments and techniques used to derive 
four SSI data sets that are employed for CCM simulations 
(Sect. 12.21 ). Finally, we discuss recent re-assessment of the 
TSI absolute value (Sect. 12.31 ) since it enters some climate 
and SSI models, and outline the perspectives of future SSI 
observations (Sect. l2~4l . 

2.1 Observation of solar irradiance 

For many centuries, the Sun has been considered as an exam- 
ple of stability and, not surprisingly, the TSI, i.e. the quantity 
of radiative power density (Wm~ 2 ) at normal incidence on 
top of the atmosphere, at a Sun-Earth distance of one astro- 
nomical unit, has been called solar constant. This term is 
still often used. Ho wever, regular space- based observations 
that started in 1978 dWillson et all Il98lb have revealed that 
the TSI varies over time scales of minutes to d ecades, and 
probably even longer (e.g. iFrohlichl 120061 12009b . 

The most noticeable variation of the TSI is a 0. 1 % mod- 
ulation in phase with the solar activity or sunspot cycle. Al- 
though quantifying such small variations is a major tech- 
nological challenge, it is strongly motivated by the desire 
to understand solar variability, and even more so by the 
importance of the TSI for terrestrial energy budget studies 



(Trenberthet al., 2009) 



'TOSCA - towards a more complete assessment of the impact 
of solar variability on the Earth's climate, http://www.tosca-cost.eu 



The TSI is the spectral integral of SSI over all wavelengths 
but its weak variability masks the fact that relative SSI varia- 
tions show a strong wavelength dependence (Fig. [TJ. In par- 
ticular, the visible and NIR bands are the least variable of 
the solar spectrum with a relative solar cycle amplitude of 
the same order as for the TSI (0. 1 %), whereas values of 1 
to 100 % are observed in the UV variations, and in excess of 
100 % in the soft X-ray range (below 10 nm). Each individ- 
ual spectral band has a markedly different impact on the ter- 
restrial atmosphere, which depends on the atmospheric pro- 
cesses affected by the given band, the amount of the spectral 
flux, and its variation. 

Figure [3] shows an overview of all SSI space experiments 
as a function of the period of observations and wavelength 
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Fig. 3. Overview of the main satellite missions that have made SSI observations at wavelengths higher than 100 nm. Short missions such as 
sounding rockets are not indicated. 



coverage (above 100 nm). Regular space-based monitor- 
ing of the solar spectrum over a broad range, covering, in 
addition to the UV, the visible and the IR up to 2.4 urn, 
started with the launch of ENVIS AT/SCIAMACHY in 2002 295 
and SORCE/SIM in 2003. Measurements in the UV below 
400 nm began several decades earlier. Most of these earlier 
measurements, however, are difficult to use for quantifying 
the solar variability, for reasons that will be detailed below. 

The longest records of SSI measurements were provided 300 
by SOLSTICE and SUSIM aboard th e UARS (Upper Atmo - 



sphere Research Satellite) spacecraft (IRottman et all 12004) . 



These instruments observed solar UV radiation between 120 
and 400nm from 1991 to 2001 and 2005, respectively. These 
measurements pointed to the impor t ance of the irradiance 305 



variat i ons in the UV (Wood s et al. l2000t IRottman et al 



Merging all UV observations into a single homogeneous 
composite record is a major challenge (DeLand et all 120041) 
that is hampered by several problems. First, the lifetime of 
most instruments does not exceed a decade. This makes 
long-term observations covering periods that exceed single 
instrument lifetimes, of prime interest for climate models, 
very difficult. A second obstacle is the differing technolo- 
gies and modes of operation of various space-based instru- 
ments. The cross-calibration of individual records is further 
hampered by the fact that overlapping observations disagree, 
and the existing data records are spectrally and temporally 
intermittent. Although missing observations can be filled in 
by using data regression based on time series of solar proxies 
such as the Mg II index, which are we l l correlated w i th UV 
variat i ons (DeLand and Cebula , 1993 ; Viereck et al. , 2001 



200 it iFlovd et all l2002l 120031) . although the solar cycle iLeanl [l997j), none of the existing solar proxies can prop 



variability of solar radiation above approximately 250 nm 
remained relatively uncertain due to insufficient long-term 
stability of the instruments (IWoods et all Il996l) . In ad- 310 



erly re produce solar irradiance in a spectral band on all time 



scales JPudok de Wit et allEooli) 



dition to SOLSTICE and SUSIM, there is the long time 
series of 200^400 nm solar UV observations by several 
NOAA SBUV instruments, w hich had underfligh t calibra- 



280 ti ons aboard the Space Shuttle (ICebula et al.l 119981) . Table 1 



of DeLand and Cebul3 ( 2012 ) summarises the measurement 315 
uncertainties for these instruments. 

Regular observations of the visible and NIR bands 
coveri ng more than o n e year started with SOHO/VIRGO- 



285 SPM (IFrohlich et all 1 19971) in 1996 at th ree selected 



bands and continued w ith th e ERS/GOME (Web er etal 



The most critical issue for all SSI instruments is the op- 
tical degradation caused by the energetic radiation in the 
space environment. Two options have been employed to 
account for instrumental degradations. The first one is 
to provide redundancy in the instrument design, by us- 
ing, e.g. a dual spectrometer setup with detectors that ex- 
perience different accumulated exposure time or by plan- 
ning redundancies in spectral channels and calibration lamps 
from which degradation corrections can be derived. This 
approach was used, for instance, for UARS/SUSIM an d 



SORCE/SIM (Brueckneretal, 1993; Harder etal.. 2005a). 



Burrows et al 



dBovensmann et all fl999 ) 
I2005a| b|). and ISS/SOLSPEC dThuillier et al 



1999), EN VIS AT/SCIAMACHY 
SORCE/SIM 



ments. 
below. 



Hard er et al 
2009|) instru- 
These SSI missions are described in more detail 



The second option is to use reference data, e.g. solar irra- 
diance measurements for long-term calibrations, as reported 
for the NOAA SBUV instruments, periodic recalibrations 
using sounding rockets, or stable external calibration tar- 
gets like selected stars, as done for UARS/SOLSTICE and 
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SORCE/SOLSTICE (iMcClintock et all 120051) . Sensitivity 
changes and degradation are strongly wavelength-dependent, 
which makes creating a properly cross-calibrated SSI record 3ao 
very difficult. One atte mpt for a composite SSI tim e series in 

330 the UV is provided bv lDeLand and Cebulal (12008b . 

Cross-calibration of different SSI records is also limited by 
the lack of realistic confidence intervals for the existing data. 
This aspect has been thoroughly investigated for the TSI and 385 
sound estimates have been obtained for time scales from days 

335 to years. No meaningful estimates, however, exist for time 
scales exceeding one decade, which are needed for climate 
studies. The situation is much worse for SSI measurements, 
whose relative uncertainties often are two to three orders of 390 
magnitude larger. 

340 In the following subsections a brief summary of newly 
available SSI observations obtained during the last decade is 
given. Also a brief description of recent developments re- 
garding the absolute values of the TSI, and its variability is 395 
provided. 

345 2.2 SSI time series 

2.2.1 SCIAMACHY and GOME 

SCIAMACHY (SCanning Imaging Absorption spectroM- 400 
eter for Atmospheric CHartographY) and GOME (Global 
Ozone Monitoring Experiment) onboard the ENVISAT and 
350 ERS-2 satellites, respectively, are atmospheri c sounders that 
measure terrestrial atmospheric t race gases (burrows et al. 



Bovensmann et al.. 1 19991: Bovensmann et all 1201 II). 405 



The primary purpose of direct solar measurements by SCIA- 
MACHY and the two GOMEs (a second GOME is flying 
on METOP-A since 2006 and a third GOME on METOP- 
B since 2012) is to Sun-normalise the backscattered light 
from the terrestrial atmosphere, which is then inverted to 4 
determine atmospheric trace gas amounts. This normalisa- 
tion does not require absolute radiometric calibration and 
cancels out degradation effects. The SCIAMACHY and 
GOME instruments have been radiometrically calibrated be- 
fore launch. 

In order to provide estimates for solar cycle variability 
from SCIAMACHY measurements (230nm-2.4um) with- 
out the need for a detailed degradation correction, a proxy 
model (hereafter referred to as SCIAMACHY proxy model) 
was developed by fitting solar proxy time series to observed 420 
SCIAMACH Y measurements ove r several 27-day solar rota- 
tion periods dPagaran et all 120091) . The pr oxies used in this 



mode l are the photometric s unspot inde x (Balm aceda et al- 
20091) and the Mg II index dWeber et all l2012b for sunspot 



darkening and facular brightening. Assuming that the fitting 425 
parameters linearly scale from solar rotations to an 1 1-yr so- 
lar cycle, one can then use the solar proxies to extrapolate 



375 beyond the lifetime of the single instrument dPagaran et al 



201 lbl) . Note, however, that this assumption might not be ac- 



nitude of the solar cycle variation. The short-term variability 
of the SCIAMACHY proxy model over several solar rota- 
tions agrees well with direct solar observatio ns from SCIA 



MACHY, SO RCE/SIM, and UARS/SUSIM dPagaran et al. 
20091 1201 lal) . Over longer periods during the descend 



ing phases of solar cycles 21 to 23, larger differences be- 
tween model and direct observations become apparent. In 
particular, SORCE/SIM data show UV changes that are 
about four times large r than the SCIAMACHY proxy model 
dPagaran et alll201 lbl Fig.EJ. 

Very recently an optical degradation model has been de- 
veloped that uses the various light paths from different com- 
binations of mirrors within SCIAMACHY. The main causes 
for the optical degradation are believed to be contaminants 
on the mirror surfaces. This new degradation model with im- 
proved calibration corrections should allow in the near future 
to derive long-term time series of SCIAMACHY measure- 
ments without the use of the SCIAMACHY proxy model. 

2.2.2 SORCE 

The Sol a r Rad iation and Climate Experiment (SORCE; 
Rottmanl 2005), launched in January 2003, has made con- 
tinuous daily measurements of SSI from 0.1 to 2400 nm 
(with missing portions of the extreme UV between 35 nm 
and 115nm), accounting for about 97% of the TSI. These 
measurements are unique compared to SSI existing data. 
TSI is also measured on SORCE by the Total Irradi- 



ance Monitor (TIM; iKopp et all 120051) . The two instru- 



ments onboard SORCE pertinent to this study are the So- 
lar Stellar Irradiance Co mparison Experiment (SOLSTICE; 
McClintock et~aill2005l) and the Spectral Irradiance Monitor 



(SIM: lHarder et alll2005 a b) 



SORCE/SOLSTICE is a grating spectrometer that mea- 
sures SSI in the UV from 115nm to 320 nm with a res- 
olution of 0. 1 nm and with an absolute calibration un- 
certainty of approximately 5 % and measurement pre- 
cision to better than 0.5 % on all time scales (Snow 
et al. 2005). SORC E/SOLSTICE is a se cond generation 



415 of UARS/SOLSTICE dRottman et all 1 1993b which acquired 



curate, and probably results in an underestimate of the mag- 430 



UV measurements from 1991 to 2001. SOLSTICE uses 
nighttime observations of stars to track and correct for 
changes in responsivity. SORCE/SOLSTICE uses two chan- 
nels to cover the spectral regions 1 15-180 nm and 170- 
320 nm. The long-term stability in the latest data version 
is about 1 % per year (M. Snow, personal communication, 
2012). 

SORCE/SIM dHarder et alll2005albl) was developed to re- 
place the longest wavelength channel (280^420 nm) in the 
original UARS/SOLSTICE and to extend wavelength cover- 
age well out into the NIR. SORCE/SIM employs a single op- 
tical element, a Fery prism, for dispersion and to focus light 
on four detectors in the focal plane. Two photodiode detec- 
tors cover the range from 200 nm to 950 nm, another covers 
the range from 895 nm to 1620nm, and an electrical sub- 
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stitution radiometer (ESR) operates over the spectral range 
from 258 nm to 2423 ran. This ESR is also used to calibrate 
the other three. Because SORCE/SIM is a prism spectrome- 
ter its resolution varies from less than 1 nm in UV to approxi- 
mately 40 nm in the NIR. SORCE/SIM began reporting daily 490 
SSI results in April 2004. It achieves an absolute calibra- 
tion uncertainty of approximately 2 % and mea surement pre- 



Merkel et al 



cision of 0.1 % or better at most wavelengths. 
d201 lb stated that the SIM long-term uncertainty in the 200- 
300 nm region is re 0.5-0.1%, from 3 10-400 nm it is re 0.2-495 
0.05%, and in the 400-1600 nm range it is better than 0.05%. 
However, they also reported the lack of independent observa- 
tions for direct validation of these estimates. Moreover, the 
values above are relative to annual changes. Therefore, they 
can be assumed valid for short-term variations, but on longer 500 
time scales, the instrument stability could be considerably 
lower. 

It is worth mentioning that the SORCE public Level 3 
data include the SOLSTICE data up to 308 nm and SIM 
data above 308 nm. However, the SIM spectra do extend 505 
down to about 200 nm. Although the SIM data between 200- 
308 nm are not publicly available yet, they have been made 
exploitable for several studies through personal communica- 
tions. The SORCE data used in our study to estimate the 
atmospheric response to SSI solar cycle variations are speci- 
fied in Table [3] They also include data received by personal 
communications with instrument team members (J. Harder, 
p ersonal comm u nicatio n, 2012). 



Harder et al 



(12009b presented multi-year SORCE/SIM 
460 trends indicating that UV variability during the declining 515 
phase of solar cycle 23 (between 2004 and 2008) was larger 
than that observed in previous cycles, and was compensated 
by trends in other bands in the visible and NIR that in- 
creased with decreasing solar activity. SORCE/SOLSTICE 
465 has also shown enhanced UV variability for the same time 
period. Solar UV variability measured by both SORCE/SIM 
and SORCE/SOLSTICE exceed the variability observed by 
UARS/SOLSTICE and UARS/SUSIM over solar cycle 22 520 
and ascending phase of cycle 23 by a factor of 3 -10 depend- 
470 ing on wavelength ( DeLand and Cebulal 2012 : Figs. 2, 4, 
and 8). 

These discrepancies with prior cycle observations and with 
SSI models have inspired new analyses and collaborations 525 
aimed at a better understanding of the potential sources of 

475 instrument degradation that might have affected SORCE in- 
struments and previous instruments as well. The studies have 
been concentrated on SSI instrument observations, capabili- 
ties, and estimated spectral irradiance uncertainties, methods 530 
of correcting for degradation, and refinement of estimated 

480 uncertainties. It has been understood that all detectors and 
optics suffer some degradation in space, largely due to expo- 
sure to solar light, and also due to hydrocarbon contamina- 
tion that dominates below 400 nm. Accordingly, new models 535 
of degradation based on total dose, rather than just exposure 

485 time, are being developed for the SORCE and other instru- 



ments. Revised data sets, which e.g. will include SIM data 
down to 240 nm, are expected out in 2013. Besides, degra- 
dation trends have also been analyzed by considering the ex- 
pected invariance of SSI over the sol ar cycle minim um. The 
latter method has been developed by IWoodsl (120121) and ap- 
plied to data during last solar cycle minimum (2008-2009) 
to estimate possible degradation trends for SORCE/SIM and 
SORCE/SOLSTICE. It consists of identifying near-identical 
solar activity levels on both sides of the minimum to derive 
corrections for instrument degradation. This analysis showed 
good agreement of the variability from moderate solar activ- 
ity level to minimum level from various measurements and 
models, from 120 nm to 300 nm for solar cycles 21 through 
24. However, as the method has about 30% uncertainty in 
variability due to the assumptions about selecting times of 
similar irradiance levels, the results may not be as accurate 
as those derived from analyses based on instrument degrada- 
tion alone. 



The analysis by IWoodsl (120 121) reduces the variability of 
the integrated UV irradiance from 200 nm to 400 nm, rela- 
tive to the measured TSI ch ange, to 110 % (Fig . I3 from the 
190% change reported by Harder et al.l (120091 ). Neverthe- 



less, to be compatible with the observed TSI changes even 
this lower amplitude of the UV variation over the solar cycle 
still requires compensation from out-of-phase trends at other 
wavelengths, in particular above re 400 nm. Other analy- 
ses of solar cycle variability suggests that the UV variability 
in the 200 nm to 400 nm range is about 60 % of the mea- 
sured TSI change dKrivova etaL 20061 : Pagaran et al. . 20091: 



Morrill et al 



201 lrT 



2.2.3 SOLSPEC 



The SOLar SPECtrum instrument (SOLSPEC; 



Thuillier et al.L 120091) is composed of three double mono- 



chromators (170-390nm, 380-850nm, 800-3000 nm) and 
a set of lamps allowing corrections for aging related to the 
harsh space environment. The SOLSPEC spectrometer flew 
several times on the Space Shuttle and its twin instrument, 
SOSP (SOlar SPectram), was placed on the EURECA 
(EUropean Retrieval Carrier) platform for 10 months (April 
1993 to January 1994). These missions have provided data 
to build the ATLAS (ATmospheric Laboratory for Applica- 
tions and Science) spectra, specifically AT LAS- 1 (March 



1992 ) and ATLAS-3 (November 1994) dThuillier et al 



2004), which are composites using UARS/SUSIM and 
UAR/SOLSTICE data from Lyman-a at 121 nm to 
200 nm, and ATLAS/SSBUV, ATLAS/SUSIM, and AT- 
LAS/SOLSPEC from 200 to 400 nm, ATLAS/SOLSPEC 
from 400 to 850 nm, and EURECA/SOSP from 800 to 
2400 nm. The ATLAS spectra were calibrated to absolute 
radiometric scale using the blackbody of the Heidelberg 
Observatory, and tungsten and deuterium lamps calibrated 
byNIST. 
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Fig. 4. Daily averages of the spectral irradiance in the 220-240 nm 
band, as measured by 6 different (and only partly overlapping) in- 
struments. Each record has been shifted vertically for easier visu- 
alisation. Extrapolations are shown in grey, observations in colour. 
One-a confidence intervals for the former are not shown here but 
are typically 0.4 x 10~ 3 (Wm^nm -1 ). 
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SOLSPEC has been up-graded for operations onboard the 
International Space Station (ISS) by implementing several 
changes in the electronics, optics, and mechanisms and by 
adding redundant components in order to generate data with 
proper degradation correction. ISS/SOLSPEC has been cal- 
ibrated to an absolute scale at the Physikalisch-Technische 
Bun desanstalt (PTB) usin g the BB3200pg blackbody radi- 



545 ator ( ISperfeld et al 



1998b . Over the whole spectral range, 
SOLSPEC accuracy is within 2 to 3 %. The SOLSPEC spec- 
trometer has been in operation since February 2008 onboard 570 
the ISS along with the SOL- ACES (SOLar Auto-Calibrating 
EUV/UV Spectrophotometers) instrument measuring below 
150 nm. When the ISS orientation allows, ISS/SOLSPEC 
records the solar spectral irradiance. Presently, data have 
been obtained during the solar minimum preceding solar 575 
cycle 24 and at some specific periods during its rising 
phase for direct comp arisons with SORCE/SIM observations 



555 dThuillier et al.L 120121) . 



2.2.4 Statistical analysis of SSI time series seo 

As mentioned above, the assessment of long-term variations 
in SSI observations can be done only by stitching together 
different records that, individually, often do not last for more 
than a few years and are almost always offset by different sss 
calibration scales. The first systematic effort toward build- 
ing such a single composit e data set for the UV was done by 



DeLand and Ceb ula (2008), who created a data set with daily 



spectra covering the wavelength range 120-400 nm for the 
time period November 1978 to August 2005. However, the 590 
instruments frequently differed in their radiometric calibra- 



Fig. 5. Modulation amplitude (top) and phase (bottom) from 
SORCE/SIM and SORCE/SOLSTICE, compared to that from other 
instruments. The SSI has been averaged over 10 nm bins. The ver- 
tical width of the shaded area reflects the dispersion of the instru- 
ments (±cj confidence interval). 



tion and in their long-term stability. One is therefore often 
left with making subjective adjustments that can dramatically 
alter the interpretation of long-term variations ( LockwoocU 
201 II) and, for example, mimic a solar cycle variation that 
does not exist. Normalisation to a s i ngle r eference spec- 
trum, as done by iDeLand and Cebu la (2008), can only par- 
tially compensate such differences in absolute calibration. 
A more objective approach consists of incorporating these 
instrumental discrepancies in the reconstruction, and then 
explicitly using them as a contribution to the overall uncer- 
tainty. A Bayesian statistical framework is ideally suited for 
this analysis. 

To illustrate this approach, we concentrate on the solar cy- 
cle (i.e. decadal) modulation, which has the advantage of 
being one of the conspicuous signatures of solar variability 
in climate records while offering sensitive means for diag- 
nosing the quality of the solar observations. 

Our analysis consists first in extending in time all SSI ob- 
servations, and subsequently in analysing their 11 -yr mod- 
ulation. First, for each instrument and for each wavelength, 
we extrapolate the SSI observations backward and forward in 
time while preserving their statistical properties with respect 
to all other observations (cross-correlation, etc.). This is 
made possible thanks to the empirical evidence for all neigh- 
boring spectral bands to evolve remarkably coherently in 
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time, on time scales of hours and beyond (iLean et all 11982 



Amblard et alll2008l) . This coherency, which is rooted in the 



strong magnetic coupling between solar atmospheric layers, 
allows us to describe all salient features of the variability in 650 
the SSI with just a few degrees of freedom. There are typ- 
ically three of them in the UV. The extrapo lation is based 
on th e expectation-maximisation technique (IDudok de Wii 



201 lb . which is routinely used to fill in missing data in cli- 
mate records. As a result, we obtain for each wavelength 
as many records as there are instruments observing it. The 
overall dispersion of the various observations is then natu- 
rally reflected by the dispersion of the reconstructions, and 655 
so no offset or trend adjustments are required. We would like 
to highlight that no proxies are used to constrain the recon- 
struction. 

The next step consists in comparing the 1 1-yr modulation 
amplitude and phase for each reconstruction and identifying 
possible changes since 1978. To date, all spectral bands and 
all solar proxies show in-phase variations with each other (up 
to a constant pass shift), regardless of the solar cycle. This is 
a strong, albeit not sufficient, indication that different spec- 
tral bands are likely to remain in phase from cycle to cy- 665 
cle. Likewise, their modulation amplitudes are unlikely to 
change. Any departure from this picture should thus receive 
special consideration. We estimate the modulation ampli- 
tudes and phases by running an 1 1-yr running mean through 
the data. Here, the phase reference is chosen to be the Mg II 670 
core-to-wi ng ratio because it is widely used as proxy for the 



for the SSI to be unusual during the last solar cycle only, we 
conclude that these observations are likely to be affected by 
instrumental drifts , in ag reement with the conclusions from 
Lean and DeLandl 020121) . This may also be the case for 
SORCE/SIM between 308 and 340 nm, although the depar- 
ture from other observations is much less significant here. 

2.3 TSI time series 

Some climate models still include solar variability only as 
derived from TSI, thereby fixing the relative spectral contri- 
bution of the solar radiation hitting the Earth. In addition, 
TSI time series also provide constraints for the empirical and 
semi-empirical models of SSI variations that are discussed 
in the following. Therefore, it seems appropriate to discuss 
here the current knowledge on TSI measurements. In fact, 
for many years, the canonical value of the average TSI was 
1365.4 ±1.3 Wm~ 2 . Now, the most accurate, and generally 
accepted, value is 136 1 ±0.5 Wm -2 ( Kopp and Leanl 2011; 



Schmutz et al., 2012). This lower value will be used in the 



data assimilation and meteorological re-analysis project like 
ERA-CLIM at ECMWF (D. Dee, ECMWF, personal com- 
munication 2012, and see http: //www.era-clim.eu/l l. 

TSI variability was alrea dy predicted i n the 1920s from 



ground-based observations (lAbbot et all 119231) . Accurate 



measurement of TSI and detection of its variability requires 
observations from space. The first report of the variable solar 



solar UV dViereck et all l200ll) . This choice, however, has irradia nce with correct amplitudes was made by|Hickey et al 



no impact on the results that follow. 

Figure |4] illustrates the extrapolation for the 220-240 nm 
band, which is im portant for ozone production (e.g. f 
Rozanov et all 2002 ). This example shows that all obser- 



1980h . Later observations differed markedly in their abso- 



vations vary in phase with the solar cycle but differ consider- 
ably in their modulation amplitude. The observations from 
SORCE exhibit a larger modulation amplitude, as already 
mentioned above. Note that the use of the Mgll index as 680 
a reference is just to fix the phase reference for comparing 
the results obtained for different cycles, without affecting the 
way the SSI is extrapolated back- and forward in time. 

Because this approach provides as many UV composites 
as there are observations, we can test whether the observa- 685 
tions from SORCE/SIM and SORCE/SOLSTICE are com- 
patible with those from other instruments, most of which 
were obtained during preceding solar cycles. The compar- 
ison is summarised in Fig. [5] which compares the modula- 
tion amplitude and phase of SORCE versus the distribution em 
obtained by the other instruments. We conclude that all in- 
struments agree remarkably well below 200 nm; at longer 
UV wavelengths, the modulation amplitude inferred from 
SORCE/SOLSTICE is systematically larger by a factor of 
two to six at all wavelengths. The simultaneous sharp drop 695 
in its phase raises doubts about the consistency between 
SOLSTICE and the other measurements, and also between 
SIM and the other data. Assuming that there is no reason 



lute values but all basically agreed in the relative amplitude 
of the TSI variations. In TableQ]we list the main space exper- 
iments that have measured TSI, together with their observed 
variabilities, which for some experiments were subsequently 
revised as further detailed below. The given numbers are bi- 
ased by the duration of the experiments and their phase rela- 
tive to the solar cycle, but the overall result is that TSI varia- 
tions are observed to be on the order of about 0.5 % standard 
deviation from the mean value. 



Lee et al 



(1995) estimated the absolute accuracy of 
the NIMBUS7/HF instru ment to be 5% and that 
of ERBS/ERBE 0.2%. IWillsorJ ( 1 19791) expected his 
SMM/ACR IM-I experiment to rema in within 0. 1 % for at 



least a year. iFrohlich and Leanl (11998b state that the absolute 
measurements of the early radiometers are uncertain to about 
0.4%, which corresponds to 5.5 Wm~ 2 . The main reason 
for this relatively large spread was due to the uncertainty in 
the aperture area. When the technology for determining aper- 
ture area improved the agreements among the measurements 
improved. However, the SORCE/TIM expe riment proved 
to be a new outlier. Lawrence et all (120031) claim an un- 
certainty of 0.5 Wm~ 2 , i.e. accurate to 350 ppm. Because 
SORCE/TIM is 4.5 and 5 WitT 2 below SOHO/VIRGO and 
ACRIM/ACRIM-III, respectively, th e uncertainties given b y 
the instrument teams do not overlap dKopp and Leanl l201ll) . 
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Table 1. Summary of main TSI measurements from space. 



Spacecraft 


Instrument 


Start 


End 


Mean and st.dev. 

(Wm- 2 ) 


Reference 


NIMBUS 


HF 


Nov 1978 


Dec 1993 


1372.1 ±0.8 


Hickev et al. f 19801 


SMM 


ACRIM-I 


Feb 1980 


Jul 1989 


1367.5 ±0.7 


Willson (19791 


ERBS 


ERBE 


Oct 1984 


Mar 2003 


1365.4 ±0.6 


Lee et al. f 19871 


UARS 


ACRIM-II 


Oct 1990 


Nov 2001 


1364.4 ±0.5 


Papetal. ( 19941 


SOHO 


VIRGO 


Feb 1996 




1365.7 ±0.6 


Frohlich et al. ( 19971 


ACRIM-sat 


ACRIM-III 


Apr 2000 




1366.2 ±0.7 


Willson (20011 


SORCE 


TIM 


Mar 2003 




1360.9 ±0.4 


Lawrence et al. (20001 


PICARD 


PREMOS 


Nov 2010 




1360.5 ±0.4 


Schmutz et al. (20091 



The PREMOS experiment on the French satellite 
PICARD, which was launched in July 2010, has con- 
tributed to the understanding of the instrument offsets, by 
confirming th e lower TS I value initially reported by the 
SORCE/TIM dKopp et all 1 20051) . The radiometers of the 745 



PICARD/PREMOS experiment have been calibrated in two 
705 different and independent way s. The first is a calibr ation in 
power response as reported by ISchmutz et al.l (120091) . In ad- 
dition, the TSI radiometer of PICARD/PREMOS is, so far, 
the first space instrument that has been calibrated in irradi-750 
ance in vacuum. This was done at the Total solar irradi- 
710 ance Radiometer Facility (TRF) located at the Laboratory 
for Atmosph eric and Space Physics (LASP) in Boulder, Col- 
orado, USA (IFehlmann et all l2012h . The irradiance calibra- 
tion is accurate to 330 ppm. PICARD/PREMOS agrees with 755 
SORCE/TIM to within 0.4 Wm- 2 with PICARD/PREMOS 
715 being lower (ISchmutz et all I2012I) . Thus, the new experi- 



ments are well within their common uncertainty range and 
the uncertainty difference between independent measure- 
ments of the TSI has decreased by a factor of ten. 760 

The new experiments and the work carried out by the in- 
ternational community, which included realization of new 
facilities, ground-based tests, and collaborations aimed at 
identifying, quantifying, and verifying the causes of the dis- 
crepancy between the TIM and older TSI instruments, have 765 
ultimately led to the understanding of the instrument off- 
sets. In particular, the characterisation of the SORCE/TIM, 
ACRIM/ACRIM III, and SOHO/VIRGO witness units at 
TRF, and the calibration of PICARD/PREMOS, led to bet- 
ter characterise instrumental effects and to r esolve the source 770 
of the discrepancy among TS I observations ( Kopp and Leant 



2011; Schmutz et all 12012 ). Instruments such as PM06 



and ACRIM type having a view-limiting aperture in front 
and a smaller precision aperture that defines the irradiance 
area have a large amount of scattered light within the in- 
strument. This additional light is not fully absorbed by the 
baffle system and produces scattered light contributing to ex- 
tra power measured by the cavity. Scattered li ght was one 



of the potential systematic errors suspected by iButler et al 
(2008). Subsequent ground testing involving the different 



roneously high values of other TSI instruments. Accord- 
ingly, new stray light corrections have recently been assigned 
to ACRIM/ ACRIM-III (based on spare instruments) and its 
stray light contribution is indeed of the order as the observed 
differences. For SOHO/VIRGO the scattered light issue was 
not the reason for its discrepant reading. VIRGO is trace- 
able to the World Radiometric Reference (WRR). New anal- 
yses, h owever, have revealed that the WRR has a systematic 



offset (IFehlmann et all 120121) . In particular, the WRR off- 



set produced approximately the same systematic shift as the 
scattering error. Thus, scattered light could be the reason of 
discrepant readings in both VIRGO and WRR instruments. 

It is worth emphasizing that the key issue in irradi- 
ance measurements is the traceability of the instrument, 
which is necessary to meet metrological requirements. The 
various instruments have undergone increasingly precise 
post-launch corrections to meet such requirements; PI- 
CARD/PREMOS was the first to achieve complete traceabil- 
ity. However, there is now a consistent evaluation of the 
TSI from four instruments: SORCE/TIM, SOHO/VIRGO 
(corrected), ACRIM/ ACRIM-III (corrected), and from PI- 
CARD/PREMOS. At the same time, TSI changes can be 
measured to much higher precision than absolute TSI val- 
ues. On short time scales, relative measurements are accurate 
to a few ppm on a daily average. On longer time scales of 
years and tens of years, the stability of the measurements are 
much more difficult to evaluate. Claims of stabilities of less 



instrument teams verified scattering as the primary cause of 
the discrepancy between the TIM measurements and the er- 7eo 



than 100 ppm over ten years (IFrohl ich 2009) are most likely 
too optimistic. A more realistic estimate comes from com- 
paring independent composites that have been constructed. 
Over the time of the last solar cycle these agree to within 
about 0.2 Wm" 2 or about 20ppmyr _1 . For pre-1996 mea- 
surements even higher uncertainties for the systematic drifts 
have to be adopted. Despite these conservative assessments, 
the TSI time record is a factor of ten more accurate than any 
SSI observation. 

2.4 Discussion of SSI observations 

While the analysis of possible degradation trends for 
SORCE, SCIAMACHY, and other missions may help to re- 
solve some of the differences in existing long-term time- 
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series of measurements, to advance our understanding of SSI 
variability new and improved, upgraded observations are re- 
quired in the coming years. To this purpose, the next gener- 835 
ation SIM instrument built for the NOAA/NASA Total So- 



2012) 



lar Irradiance Sensor (TSIS) mission (Ca halan et al. 
includes many design improvements for reducing noise and 
improving in-flight degradation tracking. The TSIS/SIM in- 
strument is currently undergoing laboratory calibrations us- 
ing radiometric technology similar to that employed to re- 
solve the source of offsets among various TSI instruments 

J 1 1 L 840 

( Kop p and Leanl 1201 II) . 

The TSIS mission might be launched in 2016. It is highly 
unlikely that TSIS/SIM and SORCE/SIM observations will 
overlap in time due to the expected lifetime of the SORCE 
batteries. Nevertheless, even without any overlap in time, S45 
the next generation measurements will help to better under- 
stand the performance of previous instruments and SSI solar 
cycle variation. In addition, several satellite missions dedi- 
cated to the observation of the Earth's atmosphere (GOME-2, 
SBUV/2, OMPS) will also provide information on part of the 

850 

SSI spectra during the possible gap in observations from the 
instruments specifically designed to SSI measurements. 



3 Models of SSI variations 

855 

Longer, uninterrupted, more stable and reliable observational 
time series are critical for understanding the origin of the 
differences between SSI measurements and improving our 
knowledge of SSI variations. However, the physics of the 
underlying processes also needs to be understood better, in 860 
order to facilitate the construction of more realistic models of 
SSI variations. Such models are particularly crucial since cli- 
mate studies including stratospheric chemistry urgently need 
long and reliable SSI data sets for realistic simulations. Ac- 
quiring sufficiently long SSI time series is a long process andses 
making them more reliable requires flying multiple new in- 
struments, which will not happen for some time. Even when 
such time series do become available, they can only be ex- 
tended into the past or future with the help of suitable mod- 
els. 870 

Although considerable progress has been ma de in 



19901 Chaomanetal.. 19961 Frohlich and LeanL 19971 


Fligge et al.. 


2000; Preminger et all l2002t lErmolli et all 
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Crouch et al. 
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modelling the SSI is more difficult, leaving considerable 
room for improvement. Here we describe recent progress in 
SSI modelling. We discuss the mechanisms responsible for 
the irradiance variations (Sect. 13. It . with special emphasis a 
on the possible differences in the spectral response of the 
irradiance to the various solar magnetic features. We then 
describe the basic principles and key components of the 
models and review five current SSI models (Sect. 13. 2t that 



are available for climate studies and that were employed 
for our estimate of the atmospheric response to SSI solar 
cycle variations. Finally, we compare the models to each 
other, confront them with the available observational data 
presented in the previous section and discuss remaining 
uncertainties and open issues (Sect. 13. 3t . 

3.1 Mechanisms of irradiance variations 

Although various mechanisms have been proposed to explain 
the variation of solar irradiance, it is now accepted that ob- 
served variations in TSI (i.e. over the last 3.5 solar cycles) 
are predominantly caused by magnetic features on the solar 
surface. We cannot rule out that on longer time scales other 
mechanisms play a significant role, but this is beyond the 
scope of this paper. 

Empirically it has been known for a long time that mag- 
netic features on the solar surface are generally either dark 
(sunspots, pores) or bright (magnetic elements forming facu- 
lae and the network) when averaged over the solar disk. Two 
questions arise from this observation: why are some flux 
tubes (the theoretical concept used to describe faculae and 
sunspots) bright, while others are dark? What happens to the 
energy flux blocked by sunspots (or equivalently, where does 
the excess energy emitted by faculae come from)? 

The strong magnetic field within both small and large 
magnetic flux tubes reduces the convective energy flux. The 
vertical radiative energy flux in the convection zone is com- 
paratively small and cannot compensate a reduction in con- 
vective flux. This leads to a cooling of magnetic features. 

The magnetic features are evacuated due to the large in- 
ternal magnetic pressure and horizontal balance of total (i.e. 
gas plus magnetic) pressure. Hence, these evacuated mag- 
netic structures are also heated by radiation flowing in from 
their dense and generally hot walls. This radiation effi- 
ciently heats features narrower than roughly 250 km, making 
them brighter than the mainly field-free part of the photo- 
sphere, especially when seen near the li mb where th e brigh t 



walls are best visible dSpruitl 119761: iKeller et all 120041) 



For larger features, the radiation does not penetrate most 
of their volume (the horizontal photon mean free path is 
roughly 50-100 km), so that features greater than roughly 
400 km in diameter remain da rk (pores and sunspots); cf. 
Grossmann-Doerth et al.l (119941) . 

What happens with the energy t hat gets blocked by 
sunspots? According to ISpruij (119821) . this energy gets re- 
distributed throughout the convection zone and is re-emitted 
again slowly over its Kelvin-Helmholtz timescale, which ex- 
ceeds the lifetime of sunspots by orders of magnitude. Sim- 
ilarly, the excess radiation coming from small flux tubes 
(which act as leaks in the solar surface, since these evacuated 
features increase the solar surface area from which radiation 
can escape) is also taken from the heat stored inside the entire 
convection zone. 
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Total irradiance is simply the integral of SSI over wave- 
length. It is in the nature of integrals that quite different 
functions of wavelength, i.e. different SSI variations, can 
lead to the same variation in TSI. The differences in the rel- 
ative change in irradiance at various wavelengths is given by 
three effects: (1) The relative sensitivity of the Planck func- 
tion to temperature increases rapidly with decreasing wave- 
length. (2) Radiation at the various wavelengths is emitted 
at different heights in the solar atmosphere. This influences 
SSI because the contrast of magnetic features relative to their 
non-magnetic surroundings is height- and hence wavelength- 
dependent. (3) At very short and at very long wavelengths 
(EUV and radio), the radiation comes from the upper transi- 
tion region and corona, where the brightest commonly found 
sources are complete loops rather than just the loop foot- 
points (the flux tubes) as at almost all wavelengths in be- 
tween. 

Point 3 refers to wavelengths we do not consider here be- 
cause they interact mainly with the uppermost regions of the 
Earth's atmosphere (mesosphere and above) and hardly con- 
tribute to TSI at all. Point 2, however, is important, since 
in general the temperature in magnetic features drops more 
slowly with height in the solar photosphere and increases 
much more rapidly with height in the chromosphere than the 
solar average. This means that radiation emitted at higher 
levels in the solar atmosphere, i.e. in the UV and in the cores 
of spectral lines, displays larger changes. Together, points 1 
and 2 favour the UV to exhibit larger variations than in the 
visible and NIR. In addition to the fact that the continuum ra- 
diation comes from greater heights and the Planck function 
shows a greater temperature sensitivity, the density of spec- 
tral lines per wavelength interval also increases very rapidly 
towards shorter wavelengths. 

Now, can magnetic features be dark in the visible, but 
bright in the UV? This is in principle possible, if the mag- 
netic feature is cool in the deep atmosphere, but hot in the 
upper photosphere. For instance, pores qualitatively show 
such a temperature profile, although there are no simultane- 
ous high-resolution observations in the visible and the UV to 
decide if the temperature gradient is sufficiently extreme to 940 
produce such an effect. But pores are relatively short lived 
and too few in number. Hence, any long-lasting global dim- 
ming in the visible at times of high activity, such as that 
shown by the SORCE/SIM measurements, can only be pro- 
duced by the small magnetic elements, specifically those in 945 
the network (these are much more numerous and more evenly 
distributed). However, such small-scale magnetic elements 
are unlikely candidates to produce a decrease in the visible 
irradiance along with increased UV irradiance. Firstly, mag- 
netic elements in the networ k are bright e ven in the contin- 950 



935 uum and at disk cent re (e.g. iKobel et al.L 1201 II) . Secondly. 



Rohrb ein et al.l d201 11) have shown that the darker than aver- 



age appearance of some magnetic elements is largely due to 
spatial smearing of the observations (although there may be 
some darkening due to the inhibition of convection around 955 
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Fig. 6. Brightness of the quiet Sun (QS; top panel) and bright- 
ness contrasts of bright active components (b ottom panel) in dif- 
ferent models : SRPM jFontenla et ail l201ll. green lines) , COSI 
JShapiro et all |2010|. purple) and SATIRE dUnruh et al.l 1 19991 , 
blue). Note that contrasts depend on the position on the disc. Shown 
are averages over the entire disc under the assumption of a homo- 
geneous spatial distribution of features. 



magnetic features). Thirdly, there are also many spectral 
lines in the visible, which brighten significantly in magnetic 
elements and counteract any darkening in the continuum. 
Finally and most importantly, magnetic elements near the 
limb are always rather bright, so that averaged over the so- 
lar disk small-scale magnetic features are expected to lead 
to a brightening. All this implies that the measurements by 
SORCE/SIM of a change in th e irradiance a t visib le wave- 



lengths out-of-phase with TSI dHarder et al.L 12009b are not 



compatible with the surface magnetic field as the source of 
the SSI variations in the visible. However, the SSI variations 
in the visible are an important contributor to TSI variations 
(Fig.|2j and the above inc ompatibility woul d be strongly con- 
tradicted by the result of iBall et al.l d2012l) that 92 % of TSI 
variations are reproduced by the evolution of the magnetic 
field at the solar surface. 
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3.2 Solar irradiance models u 

It was noticed soon after the beginning of routine monitor- 
ing of TSI from space, that changes in TSI were closely re- 
lated to the evolut ion of different brightness structures on the 
visible solar disc dFoukal and Vernazzaill979l:IWillson et al.L c 



1981 



Oster et al.l Il982t Eddy et al.lll982llFoukal and Lead. 



1986). These brightness structures (such as sunspots, pores, 
faculae, plage and network) are manifestations of the solar 
magnetic field emerging at the Sun's surface. Thus their 
evolution in a global sense (and without looking too closely 102 o 
at the details of the temporal evolution) can be relatively 
well represented by different, typically disc-integrated, prox- 
ies of solar magnetic activity, such as the sunspot num- 
ber or area, plage area, the solar radio flux at 10.7 cm 
(fl0.7), the Mg II core-to-wing index, or the Ca II K, 25 
line. This has widely been used in the oldest (but still 



widely deployed) irradiance models (e.g. Don nelly et al 



1982; Foukal and Lean, 1990; Chapman et al 
Frohlich and Leani 1 19971: IFligge etall 1 1998 



Preminger et al. 12002 ; Jain and Hasanl 20041: Pagaran et al 



1996 



2012 



Leani 



20091) . In these models, the measured irradiance variations 
are fitted via a set of activity proxies through multiple re- 
gressions. 

The success and the limitations of the regression methods 
in accounting for measured TSI variations (on time-scales of 035 
days to years) gave a strong impetus to the development of 
more sophisticated and physics-based models. Such mod- 
els consider contributions of different brightness structures 
to the irradiance change separately. Thus the solar energy 
output is the sum of the fluxes emerging from all the fea-, 4o 
tures observed on the solar visible surface (corresponding to 
the solar photosphere); the number and type of disk features 
accounted for depends on the model. Usually these models 
require two prime ingredients: (1) the surface area covered 
by each photospheric component as a function of time, ando45 
(2) the brightness of each component as a function of wave- 
length and often also of the position on the solar disc. 

(1) The surface area coverage by each photospheric com- 
ponent (i.e. sunspot umbrae and penumbrae, pores, fac- 
ulae, plage, network etc.) is derived from observations.1050 
These could be disc-integrated data (e.g. sunspot number, 
sunspot area, facular or plage area, Mg II in dex, cosmo- 



genic isotope data etc.; such as in the models by Fligg e et al 
19981 ISolanki and Unruhl 1 19981: IFligge and SolankiT 



2000; 



Krivova et alll2007ll2010tlshapiro et all 1201 ltlVieira et al 



20111) or spatially resolved observations ( e.g. in the mod 



els bylFligge et allbOOOtlErmolli et al.[l2003^Krivova et al 



Ball et al., 2011) 



2003[ IFontenla et all |2004 l201ll iFontenla and Hardei 



2005i IWenzler et al.l B004 120051 120061 lUnruh et all 12008 



(2) The brightness of individual photospheric components 
is calculated from semi-empirical models of the solar at- 
mospheric stru cture using various r adiative transfer codes, 
such as SRPM ( IFontenla et alll201 ll and references therein), 



ATLAS9 dKuruczl Il993l) in SATIRE dUnruh etall Il999t 



Fligge and Solank 



dHaberreiter et al.l 12008; Sha piro et al.l 12010b . Suchcompu 



Krivova et al., 2003) or COSI 



tations are significantly complicated by the departures from 
the local th ermodynamica l equilibrium (LTE) conditions 
(Sect. 2.6 of iRuttenl 120031) in the solar atmosphere, as well 
as by the temporal and spatial bifurcati ons of the tempera 



ture and density in the solar atmospher e (ICarlsson and Stein 



19971 lUitenbroek and Criscuolil 1201 ll) . High resolution 3-D 
models (e. g. Socas-Navarrol 1201 ll) and 3-D MHD simula- 
tions (e.g. IVoglerl 120051) of the solar atmosphere are grad- 
ually becoming available. However, current models cannot 
yet re produce available observations over the entire spec- 
trum ( lAfram et all 1201 lb . Thus, at present semi-empirical 
1-D models of the solar atmosphere are a de-facto standard 
for calculating irradiance variations. Although these models 



do not account for the spatial structure and temporal vari- 
ability of the solar atmosphere, they can be adjusted to cal- 
culate t he solar spe c trum and its variabil i ty wit h high ac- 



2000 1 curacy dUnruh et all 1 19991 IFontenla et all I201U and refer- 



ences therein). These models do not necessarily catch the 
average properties of the inhomogeneous solar atmospher e 



( lLoukitcheva et all 120041 lUitenbroek and Criscuolil 120111) . 



but, if validated and constrained by the available measure- 
ments, they can be considered as a reliable and convenient 
semi-empirical tool for modelling SSI variability. 

Presently, a wide range of solar atmospheric mo dels with 



differ e nt degrees of c o mplex it y is a v ailable (e.g. | Ku rucz. 



1993 



Fontenla et all 1 19991 120091 1201 1 
~2 



1999; Penza etal. 200 



2004T 

Avrett and Loeseil 120081 IShapiro et all 120101) . most of 



Uitenbroek, 2001 



Unruh et al. 
Kuruczl 12005 



which go back to lVernazza et al.l (119811) . 

The spectra of the quiet Sun calculated with three different 
models are shown in Fig. [6] (top panel). All plotted spectra 
are in reasonable agreement with each other, though models 
used in COSI and SATIRE are a bit closer to each other than 
to SRPM. The bottom panel of Fig. [6] shows the flux differ- 
ences (i.e. the contrasts) between bright active components 
and the quiet Sun calculated with different models. The cal- 
culations were done assuming that the active regions cover 
the entire solar disk. To calculate the variability, the con- 
trasts have to be weighted by the filling factors (surface area 
coverage), which are specific for every model. Thus, differ- 
ent magnitudes of the contrasts do not necessarily imply that 
the models produce different variability. At the same time, 
the spectral profile of the contrasts determines the depen- 
dence of the variability on the wavelength. Figure [6] shows 
that the (disc— integrated) spectral contrasts of facular, plage, 
and network regions calculated with SATIRE and COSI are 
very similar and positive over the whole spectral range. On 
the other hand, the SRPM contrasts are much lower in the 
visible compared to other models and are negative at some 
wavelengths. This is not very surprising, since the SRPM 
contrasts show n in Fig. |6] are derive d from the most recent 
versions of the IFontenla et al.l ( 1999 ) model family that were 
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3.2.1 NRLSSI 
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Fig. 7. Ratio of SSI variability to TSI variability between 200 nm 
and 1800 nm in bins of 20 nm below 400 nm, 40 nm in the 
range 400-800 nm, and 100 nm in the IR above 800 nm. Shown 
are SORCE/SOLSTICE (orange) and SORCE/SIM (red) measure- 
ments between 2004 and 2008, as well as the NRLSSI (grey)„ 110 
SATIRE-S (blue) and COSI (purple) models between the maximum 
and minimum of cycle 23. 



subsequently revised to account for some observations of the 
solar atmosphere, as further detailed below. 

11 

The advantage of employing radiances computed from the 
semi-empirical model atmospheres is that they allow compu- 
tations of solar irradiance at different wavelength (i.e. spec- 
tral irradiance). This is not straightforward with the regres- 
sion models, since in this case the regression coefficients, , 
need to be estimated from observations or alternatively, ir- 
radiance changes at individual wavelengths need to be some- 
how scaled from the TSI changes or changes in the irradiance 
at some other (known) wavelength. This latter technique is 
also partly used in the UV by the NRLSSI and SATIRE mod- 
els that are described below. 



The Naval Research Labora t ory Solar Sp ectral Irradiance0 
(NRLSSI; iLean et all Il997t iLeanl |2000|) uses the photo- 
spheric sunspot index derived from sunspot area records to 
describe the evolution of sunspots in time, and Mg II, Call 
and fl0.7 disk-integrated indices to represent facular bright- 
ening. 

Below 400 nm, the spectral irradiances are derived from 
UARS/SOLSTICE observations through a multiple regres- 
sion analysis with respect to a (SOLSTICE) reference spec- 
trum. The regression analysis includes a facular brighten- 



ing and a su nspot blocking component (see lLean et al.LI 1997 



Lean, l200d for more detail). It is known that the long-term 
stability of the UARS spectral instruments (both SOLSTICE 
and SUSIM) was not su fficient to trace solar cycle variability 
above roughly 220 nm (IWoods et al. Il996h . Therefore the 



coefficients are derived from the rotational variability so as to 
avoid any long-term instrumental effects. This approach thus 
assumes that spectral irradiance changes show the same lin- 
ear scaling with a given proxy on rotational as well as cycli- 
cal time scales. 

Above 400 nm, the facular and sunspot con- 
trasts are largely based on the contrasts presented in 
Solanki and Unruhl d 1998b - They have been scaled to ensure 
that the overall (wavelength-integrated) solar cycle change 
due to sunspots (viz faculae) agrees with the bolometric 
value for the sunspot blocking (viz f acular brightening 



e-g- 



Frohlich and Lea 



derive d from TSI modelling (see 
20041) . 

The quiet-Sun spectrum in NRLSSI is effectively a com 
posite of U ARS/SOLSTICE obs ervations (below 4 00 nm) 



SOLS PEC dThuillier et all Il998h and the model bv lKurucz 
( 1991). This composite has been scaled so that its integrated 
flux corresponds to a TSI value of 1365.5 Wm~ 2 . 

Compared to the SATIRE and COSI models (described be- 
low), the NRLSSI model shows lower variations on solar- 
cycle and longer time scales between 250 and 400 nm 
(Figs. [2] and |7). This is mainly because the regression coef- 
ficients are derived from rotational variability only. Below 
250 nm (Figs. and H) and in the visible at 400-700 nm 
(Figs. |2] and H) the NRLSSI, COSI, and SATIRE models 
agree well with each other. 

In the visible and IR range, the facular contrasts in the 
NRLSSI model re l y on t he facular model atmosphere by 
Solanki and Unruh ( 19981) . which is an earlier version of 



the model bv lUnruhet al.l (119991) currently used in SATIRE. 
Thus although the NRLSSI model does show a weak out- 
of -phase variability in the IR around 1600 nm (Fig. O, inte- 
grated over the range 1000-2430 nm (Fig. |2]i the modelled 
irradiance is in phase with the solar cycle. The SATIRE, 
COSI and OAR models all display reversed variability in this 
range. Such a reversed behaviour is also suggested by the 



"http://lasp.colorado.edu/lisird/missi/ 
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SORCE/SIM data, although the magnitude of the measured 
SIM variability is much stronger than in the models. 

Integrated over all wavelengths, the NRLSSI irradiance 
(i.e. the TSI) is in good agreement with the measurements, 
although it does not quite reproduce the co mparatively low 
TSI l evel during the last minimum in 2009 (IKopp and Lead . 
201 lb . as indicated by the PMOD com posite and also in- 
dependently found by iBall et alj d2012l) with the SATIRE 
model. This is because the Mg II index employed to describe 
the evolution of the bright component with time does not fol- 
low exactly the shape of the TSI variation over the last cycle 
(although the differences appear to be within the long-term 
measurement uncertainties; M. Snow, personal communica- 
tion, 2012). 

3.2.2 SATIRE-S 

SATIRE-S0 (Spectral And Total Irradiance REconstructions 
for the Satellite era) belongs to the class of models using 
semi-empirical model atmospheres to calculate brightnesses 
of different surface features and full-disc solar images to 
describe the surface area coverage by these components at 
a given time. The full-disc data used in SATIRE are mag- 
netograms and continuum images. So far, ground-based 
KP/NSO (Kitt Peak National Solar Observatory; 1974-2003) 
and space-based SOHO/MDI (Michelson Doppler Interfer- 
ometer; 1996-2010) and SDO/HMI (Helioseismic and Mag- 
netic Imager; since 2010) data have been employed. The fol- 
lowing model atmospheres are used: the quiet Sun model by 
Kuruca d 199 lb with an effective temperature of 5777 K, sim- 
ilar but cooler m odel atmospheres for umbra and penumbra, 
and model P of Fontenla et al. ( 1999b slightly modified by 
Unruh et al. ( 19991) to achieve better agreement with obser- 
vations in the visible and near-UV. Since the ATLAS9 code 
uses the LTE approximation, which is known to fail in the 
UV, fluxes be low 270 nm are re-sc aled using UARS/SUSIM 



observations (Krivova et al., 2006) 



The SATIRE-S modelled variability was found to be 
in very good agreement, on both rotatio nal and cyclic 



time s cales, with the PM OD TSI composite dWenzler et al 



20091; IBall et al 



(Ball et al 



2011 



2012), SORCE/TIM TSI measurements 



2012b. UARS/SUSIM s p ectral irradiance 



5- 

dKrivova et all 120061 120091: lUnruh et all 12012b . as well 
as with UARS/SUSIM and UARS/SOLSTICE Ly-a mea- 
surements (Figs. [2] and [8}. On rotational times scales, 
good agreement is also found between the SSI provided by 
SATIRE-S and the SORCE/SIM and SORCE/SOLSTICE 



measurements (Unruh etal, 2008, 2012; Balletal., 2011 



see also bottom panel of Fig. [8]l and between SATIRE- 
S and SOHO/VIRGO observations in three spectral chan- 
nels in the near-UV, visible and near-IR (IKrivova et al. , 
2003). Due to strong sensitivity trends, VIRGO spectral 
data could not be used on longer time scales. On time 
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Fig. 8. Solar UV irradiance between 220 and 240 nm calculated 
with NRLSSI (black), SATIRE-S (blue) and COSI (magenta), and 
measured with UARS/SUSIM (darker green), UARS/SOLSTICE 
(light green), SORCE/SOLSTICE (orange) and SORCE/SIM (red). 
The pale green shading in the upper panel marks the period when 
the sensitivity of the UARS/SUSIM instrument (and thus the flux) 
changed, so that a s hift w as applied to the data before that (see 
IKrivova et ail I2006I . I2009I for details). Top panel shows 3-month 
smoothed values over the period 1993-2009, while the bottom panel 
is limited to the period when SORCE was in operation, i.e. af- 
ter 2003, and shows daily values, except for the COSI model, for 
which only yearly averages are available. All quantities are nor- 
malised to the corresponding mean values during the last activity 
minimum (2009). 
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scales longer than a few months, the SATIRE-S trends, how- 
ever, diverge significantly from those shown by SORCE/SIM 
and SORCE/SOLSTICE: the change in the UV is signif- 
icantly weaker in SATIRE-S, and the inverse solar cyclers 
change in the visible is not repro duced ( Figs. [7] a nd |8j 



Ball et al.11201 UlUnruh et aUl2012l) . Interestingly, |Balletal 



(2011) have shown that when integrated over the spectral 
range of SIM (200-1630nm) and corrected for the miss- 
ing UV and IR wavelengths, SATIRE-S still reproduces overdo 
94 % of TSI fluctuations measured by SORCE/TIM. The 
trends of smoothed data agree to over 99 %. At the same 
time, wavelength— integrated SORCE/SIM data show differ- 
ent trends (only about 60% of the TIM changes are repro- 
duced over the SORCE/SIM life time), though the uncertain-1245 
ties in the integr ated SORCE/SIM data are quite large (see 



Ball et al.. 201 



LJ. 



3.2.3 SRPM 

1250 

The Solar Radiation Physical Modelling (SRPM) is a set of 
tools used to construct semi-empirical models of the solar at- 
mosphere and to derive solar radiance (or emitted intensity) 
and high-resolution irradiance spectra of the Sun. SRPM 
covers all levels of the solar atmosphere from the photo- 
sphere to the corona and takes, for most species, NLTE con- 
ditions into account. The calculations include a large num-1255 
ber of atomic levels and lines, as well as molecular lines and 
molecular photo-dissociation opacities. Seven components 
have been used in SRPM until recently corresponding to dif- 
ferent features that can be identified on the Sun at a medium- 
resolution of « 2 arcsec: quiet-Sun inter-network, quiet-Suni26o 
network lane, enhanced network, plage (that is not facula), 
facula (i.e. very bright plage), sunspot umbra and sunspot 
penumbra. The atmospheric models adopted for the calcula- 
tions o f each of these features are presented in 
2009). These models go back to the ones of 



Fontenla et al 



Fontenla et al 



1999L and references therein), whic h have been revised in 



subsequent studies of lFontenla et all d2004 l2006l l2009h to 
account for observational results and for extension of the 
atmospheric layers. The observational results include e.g. 
the negative contra st measurem e nts on facular reg i ons at i27o 
some wave l ength s (IFoukal et al. Il990t iTooka et all Il992t 
Wang et all Il998b ). Recently, further modifications were in- 
troduced into the plage models to take SORCE data into ac- 
count (see below). In addition, two new components and cor- 
responding model atmospheres, namely dark quiet-Sun inter-1275 
network and hot facula, were added but not yet integrated 



into so lar irradiance calculations presented in [Fontenla et al 



(201 1). For irradiance calculations, the distribution of the at- 
mospheric components over the solar disc at a given time is 
derived from full-disc solar filtergrams; the ones employed^ 
so far include observations from the Precision Solar Pho- 



INAF Osservatorio Astronomico di Roma and at Mauna Loa, 
as well as Meudon photographic images. 



Fontenla et alj d201 II) presented SRPM-based calculations 



of daily SSI at 1 nm resolution for the period 2000-2009 
using observations carried out with the PSP T telescope in 
Rome. The atmospheric model set used in Fontenla et al. 



d201 II) is essentially that of lFontenla et al. (2009) with a few 
modificati ons, e.g. to photosp heric layers of the plage mod- 



o photospti 

els (see iFontenla et all |201 lb . These calculations produce a 
decrease in irradiance at wavelengths below w400 nm dur- 
ing the declining phase of solar activity from 2002 to 2009, 
and a flux increase at longer wavelengths, simil arly to the be 



havio ur seen in the SORCE/SIM observations dHarder et al 
20091) . The UV variability derived from the model is, 
however, significantly weaker than what is measured by 
SORCE/SOLSTICE. Moreover, the observed changes of TSI 
on time scales l onger than the s olar ro tation are not captured 



by the model. IFontenla et al.l d201 ll) argue that a possible 



reason for this mismatch could be that the current set of the 
employed atmospheric models is still missing some compo- 
nents. 

3.2.4 COSI 

The COde for Solar IrradianceQ (COSI) calculates synthetic 
solar spectra for different components of the solar atmo- 
sphere. COSI returns the most important UV lines including 
the degree of ionisation of the elements under NLTE. The 
NLTE Opacity Distribution F unctions ( Hubenv and Lanz , 
1995b . implemented in COSI bv lHaberreiter et al.ld2008b . in- 
direct ly account for the NL TE effects in several millions of 
lines. IShapiro et al. (2010) showed that NLTE effects influ- 
ence the concentration of the negative ion of hydrogen which 
results in an approximately 10% change of the continuum 
level in the visible spectrum. The contrasts between active 
regions and the quiet Sun are also affected by the NLTE ef- 
fects. This, however, does not imply that only NLTE codes 
are capable of reproducing the visible solar spectrum as the 
NLTE effects can be imitated by a slight readjustment of 



the atmosphere structure (see, e.g. Rutten and Kostikl,ll982t 



tometric Telescopes (PSPT; ICoulter and Kuhnl 119941) at the 



Shchukina and Truiillo Buenol 1200 lb . The spectrum of the 
quiet Sun calculated with COSI is in good agreement with the 
solar spectrum measured by SOLSTICE (up to 320 nm) and 
SIM (above 320 nm) onboard the SORCE satellite during the 
2008 solar minimum, and with the SO LSPEC measurement s 
during the ATLAS 3 m ission in 1994 dThuillier et all I20TTI) . 
Shapiro et al.1 d2012ab showed that COSI can accurately re- 
produce the center-to-limb variation of solar brightness in 
the Herzberg continuum retrieved from the analysis of solar 
eclipses observed by LYRA/PROBA2. 

Solar spectral irradiance variations can be modelled from 
COSI spectra by weighting them with filling factors of sur- 
face magnetic features (e.g. derived from the magnetograms 



4 http://www.digidyna.com/Results2010/ 
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as in the SATIRE model, cf. lHaberreiter et al.l (120051) . or 
from the PSPT images as in SRPM). In the current ver- 
sion of the mo del, this has been d one employing sunspot 
an d 10 Be data dShapiro et aT , 201 1 ) and and PSPT images 



dShapiroetal.U2012d) . 



The variability so derived with COSI agrees well wiuV 
SATIRE-S and NRLSSI in the Herzberg continuum spec- 
tral range (Figs. [7] and HJ and in the visible (Figs. [2] and [7). 
In the near-IR, at 700-1000 nm, COSI shows a weak in- 
verse solar cycle variability (Figs. [2] and [7). This is proba- 



bly an artifact of the model, as Shapiro et al. ( 2010l) calcu 



lated active-region co ntrasts using the model atmospheres of 
Fontenla et al.l ( T999h . which do not distinguish between um- 



bra and penumbra. Besides, they employed a single model 
for both bright plage and plage. These assumptions affect 
the position of the inversion point; the point where the influ-1350 
ence of sunspots sta rts to outweigh the br ight components. 
To account for this, IShapiro et alJ (12012a) decreased plage 
and sunspot contrasts with respect to the quiet Sun. A more 
accurate approach is under development. 

1355 

3.2.5 OAR 

Full-di sc observations carried out with th e PSPT telescope in 
Rome dErmolli et al.l 1 19981 120071 l2010h are also employed 
in the OAR (Osservatorio Astronomico di Roma) model ofi 3 6o 
solar irradiance variations. Earlier results obtained with both 
a regression met hod and a semi-empirical model were pre- 



sented by, e.g 
d2009h and lErmolli et al 



Penza et al.l d2003l I2006I) . iDomingo etal 
d20l' 1). A new semi-empirical 



model, briefly outlined below, is currently under develop- 1365 
ment (Ermolli et al., 2012). 

Images taken between September 1997 and January 2012 
have been processed and segmented to identify different fea - 
tures on the solar disc as described by Ermolli et al.l (120101 ). 
The seven classes o f atmospheric features proposed byi37o 
Fontenla et al.l (I2009I) are considered. The brightness spec- 
trum of each feature i s computed using the semi-empirical 
atmospheric mo dels oflFontenla et al. ( 2009b through the RH 
synthesis code dUitenbroekl, I2OO2I) . Thus the set of compo- 
nents in the OAR mo del is essentially the same as used byms 
Fonten la et al. H), though t heir atmospheric structure is 
based on iFontenla et al. d2009), i.e. before the most recent 



modifications applied by IFontenla et al.l d201 II) . 

The SSI variation calculated with the OAR model in the 
visible (400-691 nm) over the period 2000-2012 shows ai38o 
trend that is opposite to that measured by SORCE/SIM and 
obtained with SRPM calculations. The variability in the 
UV (200-400 nm) estimated with the model is significantly 
weaker than that depicted by the SORCE measurements. On 
the other hand, the TSI variations calculated with the models 
agree well with the PMOD composite and SORCE/TIM mea- 
surements on both rotational and cyclic time scales (Ermolli 
et al., 2012), thus suggesting that there is no need for addi- 
tional components to reproduce the observed TSI changes. 



3.3 Discussion of SSI model results 

There has been steady progress in modelling SSI variations, 
and a number of models are now available that can be used 
as input for climate studies. However, the main uncertainty 
in the models concerns the wavelength range 220^400 nm, 
which is of particular interest for climate studies. 

Of the five models discussed in this section, specifically 
NRLSSI, SATIRE-S, COSI, SRPM and OAR, only one 
(SRPM) shows a behaviour of the UV and visible irradiance 
qualitatively resembling that of the recent SORCE/SIM mea- 
surements. However, it should be noted that the integral of 
the SSI computed with this model over the entire spectral 
range (i.e. the TSI) does not reproduce the measured cyclical 
TSI changes. None of the other four models, which are in 
closer agreement with each other, reproduces the behaviour 
of the UV and visible irradiance observed by SORCE/SIM 
(Figs. [2] [7] and [8]). These models are nevertheless all in 
good or fair agreement with earlier UARS UV observations, 
TSI meas urements and t he model ba s ed on SCIAM ACHY 
data (e.g. Krivova et all [20061 120091: IPagaran et all 
Morrill et alll201 laUShapiro et alll201 lHUnruh et al. 
Lean and DeLandi l2012t iBalletall 120121) . These mod 



2009; 



2012; 



els a lso agree with SORCE data on rotational times scales 
(e.g. lUnruh et all 120081: [Ball et all I20TTI iLean and DeLand , 



20121 see also bottom panel of Fig. |8). Note that good 



agreement on rotat i onal time scales was also found by 
iDeLand and Cebula d2012l) between SORCE and other spec- 
tral observations. 

When comparing various models and data to each other, 
however, one important issue often remains unacknowl- 
edged, namely the true uncertainties in the measurements, 
in particular on time scales of years and longer. As an exam- 
ple, we compare measured and modelled variability at 220- 
240 nm (Fig.©. Over the period 2004-2008, the three mod- 
els shown in this figure (NRLSSI, SATIRE-S and COSI) sug- 
gest a decrease in the integrated 220-240 nm flux of about 
1 %. The fluxes measured by the SORCE/SOLSTICE and 
SORCE/SIM instruments in this range decreased over the 
same period by roughly 4 % and 7 %, respectively. In the 
220-240 nm region, the accuracy of S ORCE/SOLSTICE i s 
considered to be better than that of SIM (Ha rder et aill2010l) . 



Currently, SOLSTICE long-term stabil ity in this spectra l 
range is estimated to be about 1 %yr~ 1 dSnow et all 120051) . 



Thus the observed trend for SORCE/SOLSTICE (4 % over 
5 yr) and its difference with the models (3 % over 5 yr), in 
fact, lies within the long-term instrumental uncertainty. The 
difference between the trends for SORCE/SIM and the mod- 
els (about 6 %) is just outside the value of 5 % over 5 yr, 
by assuming that the SIM long-term stability is compara- 
ble to t hat of SOLSTICE Although the SIM uncertainty 
stated in lMerkel et al.l (1201 ll) is lower, this assumption seems 
to be a reasonable working hypothesis for our study, be- 
cause it is known that at 220-240 SIM has poorer precision 
than SORCE/SOLSTICE (J. Harder, personal communica- 
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tion, 2012). It is worth mentioning that new assessments of 
the SORCE instrument degradation might imply a revision 
of results derived from several studies, including this one, 
which however is mainly intended to discuss the data andi445 
models available to date for the evaluation of the atmospheric 
response to SSI variations. 

Despite showing much lower variability than SORCE/SIM 
in the UV, the models also display considerable differences 
between each other in the range 250^400 nm (up to a factor 
of three, e.g. between NRLSSI and COSI; Figs. E and E}. 
As a consequence of the low response of UARS/SOLSTICB«° 
to long-term variability above 300 nm and the use of rota- 
tional variability to estimate the regression coefficients, it is 
likely that NRLSSI underestimates the changes in this range 
(J. Lean, personal communication, 2012), and can thus be 
considered as the lower limit. All other models rely on semi-1155 
empirical model atmospheres in this range, which also need 
further tests at all wavelengths. The differences between the 
models in the IR are mainly due to the lack of reliable obser- 
vations of contrasts of different (solar) atmospheric compo- 
nents and related uncertainties in the corresponding model at-i«o 
mospheres. In the near-IR (700-1000 nm), all models qual- 
itatively agree, except COSI. However, as discussed earlier, 
the inverse variability shown by COSI in this range is be- 
lieved to be an artifact of the atmospheric components and 
models adopted for the irradiance calculations. i«s 

The variability of the spectral irradiance on time scales 
longer than the solar cycle is beyond the scope of this 
paper. We note, however, that uncertainties in the long- 
term SSI reconstructions are essentially the same as the 
ones discussed in this section convolved with the uncer-i«o 
tainty in t he magnitude of the secula r change in irradiance 



(see, e.g. |_K rivova and Solankil l2013t I Schmidt et all l2012t 



Solanki et al.Ll2013L and references therein). 



4 Climate impact of SSI measurements 



The SSI data derived from recent solar observations and 
model calculations presented in the previous sections are em- 
ployed here to discuss the impact of SSI variations on theiwo 
Earth's atmosphere and climate. In particular, we present the 
motivation for simulations of the response of the Earth's at- 
mosphere to variations induced by SSI changes (Sect. I4.lt . 
We then discuss the solar induced differences in atmo- 
spheric heating rates, ozone variability, temperature, and485 
atmospheric circulation (Sect. I4.2b . These differences are 
investigated in recent simulations using atmospheric mod- 
els with standard solar forcing from the NRL SSI model 
(e.g. lSPARC-CCMVaiEoiollTaylor et afll2012l for CMIP5 
simulations). Additionally a comparison to recent atmo-1490 
spheric model simul ations using the SORCE measurements 
as solar forcing ( e.g. |Haigh et all l2010t iMerkel et al.l |201 1[ 



SORCE measurements represent the lower and the upper 
boundaries of current understanding of SSI solar-cycle vari- 
ability. Here we provide an estimate of the uncertainty in the 
magnitude of solar induced atmospheric changes in CCMs 
by comparing the impact from the two extremes in the cur- 
rent evaluation of SSI variability (Sect. 14. 3t . 

4.1 Climate modelling 

Today, the most advanced tools available for climate simula- 
tions are 3-dimensional General Circulation Models (GCMs) 
that numerically simulate the general circulation of the atmo- 
sphere and/or the ocean based on well-established physical 
principles. Most climate models that are utilised for future 



climate predictions in the 4th IPCC report (ISolomon et al 



120071) are coupled atmosphere-ocean models reaching up to 
the middle stratosphere (32 km), whereas CCMs that are 
used for future pred ictions of the stratospheric ozone layer 
in the WMO report (iWMOl 1201 II) include interactive strato- 
spheric chemistry and reach up to the lower mesosphere or 
above (80 km). The main external driving force for all cli- 
mate models is the incoming solar flux at the top of the atmo- 
sphere (TOA). Currently there are two mechanisms for solar 
radiati on influence on climate, the so-called "top-down" UV 
effect (Kodera and Kuroda, 2002) and the "bottom-up" TSI 



effect (Ivan Loon et al.. 2007). 



As discussed in the previous sections, even though TSI 
varies only by about 0.1% over the solar cycle, larger vari- 
ations of several percent occur in the UV part of the spec- 
trum, including the ozone absorption bands between 200 and 
400 nm that are responsible for the S W heating of the strato- 
sphere and ar e important for photochemical processes (e.g. 
Haighl 119941) . The cyclic variation of the incoming solar 
irradiance at short wavelengths leads to statistically signifi- 
cant ozone, t emperature and zonal wind solar signa ls in the 



1475 stratosphere ([Austin et all l2008t iGrav et al.l I2010I) . These 



Oberlander et al.ll2012l) is provided. As shown and discussed 



in the previous sections, the NRLSSI model output and the 



solar induced circulation changes in the stratosphere can in- 
duce noticeable decadal climat e changes in the lower atmo- 
sphere and at the surface (e.g. Haighl 19991 iKoderal 2002 ; 
iMatthes et all 120061: lliieson et all 1201 1[ iMatthesl |201 ll) . In 
order to account for this so-called "top-down" stratospheric 
UV mechanism, the radiation code in climate models has 
to account for spectrally resolved irradiance changes. The 
first climate models focused on tropospheric climate, thus 
solar solar forcing was represented by TSI only. As strato- 
spheric changes played a minor role for climate predictions 
in the past, climate models did not take into account strato- 
spheric processes and in particular ozone changes due to so- 
lar UV absorption. Hence, most of the SW radiation codes 
developed for use in GCMs did not consider solar irradi- 
ance for wavelengths shorter than 250 nm and employed pa- 
rameterisation using TSI as input. Solar fluxes and heating 
rates were subsequently calculated in one or two SW ab- 
sorption bands from the top of the atmosphere to the sur- 
face. In contrast, middle atmosphere models use SW radia- 



I. Ermolli et al.: Spectral irradiance and climate 



19 



tion codes specifically designed for simulations of the upper 
atmosphere. They, therefore, cover a broader spectral range 
and include more than two spectral bands in the UV and vis- 
ible. 

Studies on the performance of SW radiation codes with 
different spectral resolution showed that the observed so- 
lar temperature signal in the stratosphere can only be repro- 
duced in models that allow for the effects of spectral varia- 
tions between solar minim um and maximum dEgorova et al 



Forster et al. 



20041: Nissen et al. , 2007 ). In a recent paper, 
d201 II) examined in detail the sensitivity of a number of CCM 
SW radiation codes to changes in solar irradiance and ozone 
as well as the ability of the models to reproduce the 11- 
yr rad iative solar signal using the NRLSSI data. In their 
study dSPARC-CCMVal 120101) . the strongest solar temper- 
ature signal was found to be in the tropical upper strato- 
sphere/lower mesosphere, indicating that the direct mecha- 
nism of heating by absorption of enhanced UV radiation at 
solar maximum is well captured by the models that employ 
spectrally resolved SW radiation schemes. Models that do 
not account for SSI variations and only consider changes in 
spectrally integrated TSI cannot properly sim ulate solar in 



duced variations in stratospheric temperature (IForster et al 
201 lh. 



Today stratospheric processes are gaining a lot of inter- 
est due to their importance for climate. Not only the ef- 
fect of ozone recovery and its relationship to climate but 
also stratosphere-troposphere dynamical coupling and its 
role for predictability from days to decades have been recog- 
nised as important issues for future c limate studies (e.g. 
Baldwin et alll2007tlGerber et al.Ll2010h . Therefore, abetter 1550 
representation of the stratosphere including improved repre- 
sentation of SW heating processes as well as dynamical cou- 
pling with the troposphere in global climate models is criti- 
cally important. 

1555 

4.2 Impact of SSI variability in climate models 

The uncertainty of SSI variations in recent observations and 
models has significant influence on simulations of the cli- 
mate system, since the response of the atmosphere strongly 1560 
depends on the spectral distribution of the solar irradiance. 
The effects for middle atmosphere heating, ozone chemistry 
and middle atmospheric temperatures are examined in the 
following. 

1565 

4.2.1 Effects on atmospheric heating and ozone chem- 
istry 

As solar radiation is the primary source of energy that drives 
atmospheric as well as oceanic circulation, accurate repre-1570 
sentation of solar irradiance is of paramount importance for 
the simulation of the atmospheric temperature, composition 
and dynamics in climate models. The variability of the so- 
lar spectrum in time, and in particular over the solar cycle is 
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Fig. 9. Updated from lOberlander eTail d2012l) (a) Difference in 
zonal mean solar heating rates between 2004 and 2007 in Kday -1 
from the SIM data, and (b) difference in the global mean 2004-2007 
solar heating rate signal between SORCE and NRLSSI data. Solid 
lines denote the UV spectral range, and dashed lines denote the VIS 
spectral range. Red dashed lines in panels (a) and (b) denote the 
same differences in zo nal and global mean, as calculated and pre- 
sented in Figs. 3b, c in Oberliind er et al.l d 20 1 2h : black dashed lines 
indicate the enhanced solar signal in the VIS range, when extend- 
ing the spectral resolution of the FUBRad shortwave radiation code 
from one to six bands in the Chappuis bands. The simulations were 
performed for perpetual lanuary conditions. 



necessary for the assessment of solar influence on climate. 
The amplitude of the simulated solar signal depends on the 
spectral solar fluxes prescribed at the TOA. Differences in the 
TOA solar irradiance spectrum result in large changes in the 
heating rates calculated by S W radiation schemes or radiative 
transfer models, as has been shown by IZhong et al.l (120081) 
with differenc es of up to sa 1 . 1 K day - 1 in mid-latitude sum- 
mer. Recently. lOberlander et al.l (120 1 21) examined the impact 
of a number of different estimates of prescribed TOA solar 
fluxes on the solar response in a GCM which inc ludes a radia 



tion s cheme with enhanced spectral resolution dNissen et al 



2007} and is therefore able to accurately represent the solar 
signal induced changes. They used the NRLSSI, the SATIRE 
and the SCIAMACHY solar flux input data sets, and com- 
pared their effects on SW heating rates over the 11 -yr so- 
lar cycle using offline calculations. They also calculated the 
corresponding temperature response from perpetual January 
GCM simulations with prescribed ozone concentrations. 

The comparison revealed clear differences in SW heating 
rates for the solar minimum of cycle 22 (September 1986). 
The simulations forced with the NRLSSI reconstructions 
show the smallest solar heating rates. The use of the SATIRE 
reconstructions leads to stronger solar heating of up to 5 % 
in the middle and upper stratosphere. The SCIAMACHY 
observations slightly enhance the solar heating in the meso- 
sphere, with differences arising from the stronger solar fluxes 
in the Huggins bands of the SATIRE model and enhanced 
fluxes in the Hartley bands of the SCIAMACHY data set. 
Using SORCE (SIM and SOLSTICE for wavelengths below 
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210 nm) measurements over the period May 2004 to Novem- 
ber 2007 reveals larger changes in solar heating rates and the 
resulting temperatures in comparison to the NRLSSI data.1630 
The lower irradiance in the visible range at higher solar ac- 
tivity than at minimum activity in the SORCE/SIM data does 
not lead to a decrease in total radiative heating. 

The spectral resolution of the SW radiation scheme of 
EMAC-FUB has recently been extended towards a more de-1635 
tailed representat i on of the Chappuis bands, and an update of 
Oberlander et al.l (120121) . performed with the extended spec- 
tral resolution for solar minimum conditions in November 
2007, separated for the UV and VIS spectral changes is pre- 
sented in Fig. [9] The enhanced UV irradiance in 2004 (com-wo 
pared to 2007) in the SORCE data leads to higher SW heat- 
ing rates while lower heating rates are simulated in the visible 
spectral range (Fig. |9^). The SW heating rate change from 
2004 to 2007 is stronger when the model is forced with the 
SORCE than with the NRLSSI data by O.lSKday" 1 in th©64 5 
global mean (Fig.[2b). An enhanced sensitivity to changes in 
the Chappuis bands is found, a result which illustrates better 
the heating rate changes in the UV, visible and NIR spectral 
regions from using the SORCE data. 

The different magnitude and spectral composition of theeso 
SSI changes leads not only to a substantial alteration of heat- 
ing rates considered above but also affects photolysis rates 
which drive atmospheric chemistry and regulate the atmo- 
spheric ozone distribu t ion in the middle atmosphere (e.g. 
Brasseur and Solomon! 2005 ). The global ozone abundance^ 
is maintained by ozone production, destruction and trans- 
port by air motions. However, in the tropical stratosphere 
above ~ 30 km the ozone concentration depends primarily 
on photochemical processes, with oxygen photolysis play- 
ing a crucial role in atmospheric chemistry. Therefore, theieeo 
spectral composition and the magnitude of the SSI changes 
are of critical importance, since they define not only the 
mag nitude but also the si gn of the direct ozone response 



(Chemistry-Climate Model Validation) and the CMIP5 
simul ati ons for the upcomin g IP CC report (ITaylor et aL 
2012b . ICahalan et all d2010l) and iHaighetalJ d2010h pub 



(e.g. iRozanov et all 120021) . The solar induced net effect of 
ozone in the stratosphere depends on the competition be-^as 
tween ozone production due to oxygen photolysis in the 
Herzberg continuum (1 85-242 nm) and ozone destruction 
caused by the ozone photolysis in the Hartley band (between 
200 and 300 nm in the UV). 

1670 

4.2.2 Effects on ozone and temperature from atmo- 
spheric model simulations 

The effects of different SSI estimates on stratospheric ozone 
as deduced by CCMs or 2-D radiative-photochemical models 
simulations have been recently reported in a number of 
papers which are discussed and compared to each other 
in the following. In particular, a number of analyses 
focused on the differences in atmospheric response between 
simulations using the SORCE and the NRLSSI modeleso 
results. The NRLSSI is the data set widely used for cli- 
mate modelling purposes, including the SPARC CCMVal 



lished first the important implications of the SORCE data 
for middle atmosphere heating and therefore temperatures. 
Using simulations from a 2-D radiative-photochemical 
model IHaigh et al. I d2010h also presented effects on ozone. 



models followed ( 


Merkeletal.1 12011 


; Ineson et al 


. 2011 


Matthes . 201 1 ; Oberlander et all 12012 


; Swartz et al. 


2012). 



To better understand the atmospheric sensitivity to the range 
of SSI estimates discussed above, these recently published 
model experiments using the NRLSSI and the SORCE data 
are compared to each other in Figs. [lO]and[TT|with respect to 
the solar signal in SW heating rates, temperatures and ozone. 
Additionally, results from three CMIP5 model experiments 
using the NRLSSI data set are shown. Table |2] presents 
an overview of the models and their experimental designs. 
Details on the spectral solar forcing used in each simulation 
are given in Table [3] The purpose is to provide the reader 
with an initial comparison of the impacts of two SSI data 
sets which represent the lower and upper boundary of SSI 
variations. One has to keep in mind, though, that all models 
used slightly different experimental setups (including 
slightly different spectral ranges for the SORCE data) and 
therefore an exact comparison between them as well as 
to observations awaits common coordinated experiments. 
Moreover, we note that the solar spectral forcing from the 
SORCE satellite used by all the models is derived from a 
relatively short period, and longer model simulations with 
full solar cycle forcing would be needed in order to better 
evaluate the differences in the response. 

Response in shortwave heating rates 

Figure [10] shows the differences in SW heating rates (top) 
and temperatures (bottom panel) for the month of January in 
the tropical region, i.e. averaged between 25° S and 25° N. 
Differences are calculated between January of the year 2004 
(during the declining phase of solar cycle 23), and January 
of 2007, close to solar minimum in December 2008, us- 
ing the NRLSSI data (solid lines) and the SORCE/SIM and 
SORCE/SOLSTICE data (dashed lines). 

Most of the models listed in the upper part of Tables |2] 
and [3] performed equilibrium simulations for the years 2004 
and 2007. EM AC performed perpetual January runs only, 
while SOCOL performed transient simulations from 2004 to 
2009, with the differences shown here calculated between 
January of the year 2004 and January of 2007. Therefore 
most of the displayed differences do not provide one full 
solar cycle signal. The HadGEM3 model simulation used 
the SORCE/SIM measurements and scaled the data from 



2004 to 2007 to obtain a full solar cycle signal dlneson et al 
201 II) . GEOSCCM simulations were performed for solar 



maximum and minimum conditions as well, inferred by 
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Table 2. Overview of atmospheric chemistry-climate models and the respective experimental designs using the NRLSSI and/or the SORCE 
data to study the atmospheric impact. 



Horizontal res- Top level 
olution 



Number of Ozone interac- QBO 
layers tively calculated 



Length of simula- Equilibrium Reference 
tion simulation 



SORCE forcing comparison 



EMAC-CCM with 


T42 


O.UlhPa 


39 


no 




no 


50 yr 




peipetual Jan 


FUBRad 




















GEOS CCM (using 


2.5° x 2.0° 


0.01 hPa 


72 


yes 




no 


25 yr 




yes 


GEOS5) 




















HadGEM3 


1.875° x 1.25° 


85 km 


85 


no 




yes 


20 x 4 yr ens 




yes 


IC2-D 


9.5° lat 


0.26 hPa for 


17 for 


chem- yes 




no 


670 days 




yes 






chemistry 


islrv 














SOCOL v2.0 


T30 


0.01 hPa 


39 


yes 




nudged 


62 months 




transient run 




















(01.2004- 




















02.2009) 


WACCM 


2.5° x 1.89° 


0.000006 hPa 


66 


yes 




yes 


25 years 




yes 










CMIP5 simulations 








HadGEM2 


1.875° x 1.25° 


85 Km 


60 


prescribed 


solar 


yes 


3 historical 


ens. 


no 










induced 


ozone 




(1850-2005) 














variation 












MP1-ESM-LR (atmo- 


T63 


0.01 hPa 


47 


prescribed 


solar 


no 


3 historical 


ens. 


no 


spheric component: 








induced 


ozone 




(1850-2005) 






ECHAM6) 








variation 












WACCM* CESM1 


2.5° x 1.89° 


0.000006 hPa 


66 


yes 




nudged 


3 historical 


ens. 


no 



Ineson et al] feOl lh 
HaighetaljlWd) 



Schraner et al 




ShaDiro et al.l 


2012b) 



iMerkeletallfcOllI) 
Hardiman et al. (2012) 



H. Schmidt (personal commu- 
nication, 2012) 

iMarshet aljb012l) 



(WACCM) 



(1960-2005) 



Table 3. Overview of the spectral range of forcing for the various atmospheric chemistry-climate models employed to study the atmospheric 
impact. 



Model 



Spectral range of forcing 



EMAC-CCM with FUBRad 
GEOSCCM (using GEOS5) 

HadGEM3 

IC2-D 
SOCOL v2.0 

WACCM 



SORCE forcing comparison 

1) NRLSSI 

2) SIM (SOLSTICE <210nm) 

1) NRLSSI 

2) SOLSTICE <244 nm and SIM 244 to 1600 nm 

Solar max and solar min conditions inferred by linear fitting approach based on 
the available SORCE time series. 
SIM (only 200-320 nm) 

Extrapolated in time to represent the full solar-cycle amplitude, The estimated 
4% change is evenly distributed in the UV band. No changes in other bands. 

1) NRLSSI 

2) SIM (SOLSTICE <200nm) 

1) NRLSSI 

2) SOCOL SIM (SOLSTICE <210 nm and SIM 210-750 nm) 

3) SOCOL SOLSTICE (SOLSTICE < 290 nm and SIM 290-750 nm) 

1) NRLSSI 

2) SOLSTICE (1 16-240 nm), SIM (240-2000 nm) and SRPM (2000- 
UOOOnm) 



HadGEM2 

MPI-ESM-LR (atmospheric component: ECHAM6) 
WACCM* CESM1 (WACCM) 



CMIP5 simulations 

NRLSSI (CMIP recommended) 
NRLSSI (CMIP recommended) 
NRLSSI (CMIP recommended) 



a linear fitt ing approach based on the available SORCE 



time series (ISwartz et al.l 120121) . The differences presented 



for HadGEM3 and GEOSCCM are slightly larger than 
for the other models (and closer to each other), because 
the simulations correspond to the full solar cycle signal.1695 
Please note that SOCOL simulations were performed with 
the NRLSSI as well as co mbinations of S ORCE/ SIM and 
SORCE/SOLSTICE data dShapiro et al.l l2012bh . The 



simulations performed with three CMIP5 models listed 
in the lower part of Tables [2] and [3] represent ensemble 
mea n time series experim ents with varying solar cycle 
(e.g. lHardiman et all I2012I) . In order to be comparable to 
the other equilibrium simulations, the differences shown 
represent the solar signal from a multiple linear regression 
analysis using the monthly mean TSI or Hartley band 
irradiance as a solar proxy. Afterwards the resulting solar 
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Fig. 10. Differences for January between 2004 and 2007 averaged over the tropics from 25° S to 25° N for different atmospheric models 
using the NRLSSI data (solid lines) and the SORCE/SIM and SORCE/SOLSTICE data (dashed lines) for (a) the shortwave heating rate in 
Kelvin per day (Kday -1 ) and (b) the temperature in Kelvin (K). For details please see text. 
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signal has been scaled to the 2004 to 2007 value, i.e. 24% 
of the full solar cycle. Note that not all the models were 
able to provide SW heating rates and/or used interactively 1735 
calculated ozone. Therefore only a model subset is displayed 
in the respective figures. In addition to the lines for the single 
model realisations a multi model mean (MMM) for three 
models, namely WACCM, GEOSCCM, and HadGEM3 is 
included in Figs. [10] and QT| As expected, the largest SW1740 
heating rate response occurs near the stratopause (^ 1 hPa) 
showing the direct impact of solar irradiance variations from 
2004 to 2007. Although the strength of the response differs 
between the different models, all models show a consistently 
stronger SW heating (by a factor of up to two to three) usingi/45 
the SORCE data as compared to the NRLSSI data. For both 
NRLSSI and SORCE simulations, the GEOSCCM shows 
the strongest response with SW heating rate differences of 



up to 0.12 Kday" 1 for NRLSSI compared to 0.42 Relay" 1 
for SORCE closer to the HadGEM3 response, as both 
models are depicting full solar cycle conditions. In contrast 
to the other models, HadGEM3, the EMAC-FUB did not 
take into account solar induced ozone variations. EMAC 
therefore shows the weakest SW response in both experi- 
ments, followed by the WACCM model. The simulations 
performed with the CMIP5 models fit nicely to the equilib- 
rium simulations with the same (NRLSSI) forcing. Given 
the fact that the NRLSSI data provide a lower boundary 
for SSI variations with the solar cycle, the SW heating rate 
effects are likely enhanced in nature as compared to recent 
chemistry-climate model simulations. The SW heating rate 
differences described above may be attributed to the differ- 
ence in the radiation codes of the models, as discussed in the 
previous section, but are also due to the different model se- 
tups (for example, in simulations performed with or without 
interactive ozone, the signal will be enhanced or diminished). 

Response in temperatures 

Corresponding to the SW heating rate differences, the 
tropical temperature differences for January are largest near 
the stratopause; using SORCE, the temperature differences 
between 2004 and 2007 are three to four times larger than 
with NRLSSI. Again, the GEOSCCM and the HadGEM 
models produce the strongest response with about 2K 
and the WACCM and the FUB-EMAC model the weakest 
response. The response in SOCOL differs slightly from 
that of the other models by showing the maximum effect 
around 0.3 hPa, i.e. at a higher altitude than the other 
model simulations. In addition, the temperature response in 
SOCOL is always weaker using SORCE/SIM or NRLSSI 
data. When using SORCE/SOLSTICE spectral forcing the 
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temperature response is approximately half of the response 
obtained using SORCE/SIM, indicating a large sensitivity 805 
of the model to the SSI data sets depicting their differences. 
Although none of the other models has performed a similar 
set of simulations, it is expected that this sensitivity is 
typical for all models shown here. Please note that the 
negative temperature response in SOCOL in the lower andsio 
middle stratosphere is not statistically significant. The 
temperature responses in the CMIP5 models closely agree 
with the NRLSSI equilibrium simulations. Comparison 
of the temperature r esponse to observations over a full 
11-yr solar cycle (e.g. Austin et al. . 2008 ; Remsbergl 2008)isi5 
reveals that the signals derived with the use the NRLSSI data 
generally shows a response peaking around the stratopause, 
with altitudes of maximum response varying between 
models according to their model setups as well as the use 
of solar induced ozone variations. As presented in Chapters 
8 of the SPARC-CCMVal Report (2011), where the solar 
cycle effect on stratospheric temperature from a number of 
CCM simulations is discussed, the response to the NRLSSI 
forcing over the full solar cycle is about 0.8K, similar to 
observations. The response to the SORCE data forcing isis25 
larger than in the observations, but does not scale linearly 
with larger forcing. It was found to be higher by almost a 
factor of two in the case of the GEOSCCM simulation which 



2012) 



corresponds to a full solar cycle (ISwartz et al 
Response in ozone 



Figure QT| displays the tropical ozone response in percent 
for month of January for those models that calculate ozone 
interactively. In the lower stratosphere, ozone changes 835 
induced by the forcing with SORCE data are larger and 
remain positive as compared to those calculated with 
NRLSSI data. In the middle to upper stratosphere ozone 
differences induced by the NRLSSI data are smaller but 
still positive whereas the ozone response using the SORCEs4o 
data is negative. The height of the change from positive to 
negative ozone response varies from 5 hPa to 2 hPa between 
the models. Observations from 8 yr of SABER data ind icate 
a transition altitude of about 1 hPa dMerkel et all 1201 lb that 
is higher than for all models. This change of sign in thei845 
ozone response using SORCE data, meaning lower ozone 
during solar maximum than during solar minimum in the 
upper stratosphere, is statistically significant in all models. 
Swartz et al.1 d2012l) . who also examined the total ozone 



response and contributions from heating and photolysis,is5o 
show good agreement betwe en observations and simulations 
with the GEOSCCM model. lHaigh et alJ d2010h showed that 
the ozone decrease in the upper stratosphere and mesosphere 
is related to photochemical processes. Decreased O3 at solar 
maximum is consistent with increased HO x and O and alsomss 
leads to a self-healing effect with more UV radiation reach- 
ing lower levels, enhancing O2 photolysis and therefore 
increased O3. Additional observational evidence is needed 



to confirm the sign reversal in upper atmospheric ozone 
response. 

Dynamical changes and impact on the troposphere 

So far the discussion focused on the response of stratospheric 
heating, temperatures and ozone to the forcing with NRLSSI 
and SORCE data. The direct solar response in the upper 
stratosphere induces ind irect circulatio n changes through- 
out the stratosphere (e.g. Kurodal 2002b and also affects the 
troposphere and the surface ( Haighi 1999 ; Kodera , 2002 ; 
Matthes et al. . 20061) . The impact on the tropo s phere using 
SORCE data has been reported bv llneson et al.l (1201 II) . who 
provide intriguing but provisional results on northern hemi- 
spheric winter circulation. Figure Q~2]displays January multi- 
model mean signals for three climate models (GEOSCCM, 
HadGEM3, and WACCM) in zonal mean zonal wind as well 
as 500 hPa geopotential height using SORCE data as ex- 
ternal forcing. A stronger and statistically significant polar 
night jet during solar maximum years dominates the strato- 
sphere and reaches down to the troposphere, whereas neg- 
ative wind anomalies that are statistically significant in 
the lower stratosphere and upper troposphere only domi- 
nate equatorwards (Fig. Q~2] top panel). The zonal wind sig- 
nals correspond to a positive Arctic oscillation (AO) pattern 
with a statistically significant stronger polar vortex and a 
partly statistically significant enhanced geopotential heights 
in midlatitudes (Fig. \l2\ bottom panel) in agreement with 
observations (e.g. iKoderai 120021) . The stronger polar vor- 
tex in the stratosphere leads to a stronger positive phase 
of AO and North Atlantic Oscillation (NAO), which means 
a stronger Iceland low and higher pressure of the Azores 
and hence an amplified storm track over Europe. This in 
turn leads to mild conditions over northern Europe and the 
eastern US and dry conditions in the Mediterranean. How- 
ever, the tropospheric AO and NAO already respond to lower 
U V irradiance variations over the 1 1 -yr cycle, as is shown 



Langematz et al.l (120 131) who obtained a similar AO/NAO 



response from a transient EMAC-FUB simulation for the 
period 1960-2005 usi ng NRLSSI d a ta as T OA input. The 
signal is seen also in Matthes et al. I d2006l) . but in equilib- 
rium simulat i ons w ith NRLSSI as TOA input. As shown by 
llneson et al.l (1201 II) these patterns imply that the solar cycle 
effect on the AO/NAO contributes to a substantial fraction 
of the typical year-to-year variations and provides therefore 
a potentially useful source of improved decadal climate pre- 
dictability for the Northern Hemisphere. Note that the re- 
sponse is regional and is negligable on the global average. 
However, a caveat is that the 1 1 -yr solar cycle variability can- 
not be forecasted into the future (on a daily, yearly or decadal 
time scale). 

Even though the solar variability on time scales longer 
than the 1 1-yr solar cycle is beyond the scope of this paper, 
we should note here again that the Sun is the fundamental 
energy source of the climate system. As such, the low so- 
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Fig. 12. Multi Model Mean (MMM) of the HadGEM3, the GEOS CCM, and the WACCM model for (a) the zonal mean zonal wind 
differences between 2004 and 2007 in January using the SORCE forcing, and (b) the corresponding geopotential height differences at 
500 hPa. Coloured areas provide 95% statistically significant signals. 



lar activity in the past few years (compared to the previous 
6 solar cycles), and its possible implications for future cli- 
mate evolution has attracted the attention of both scientists 



and the public (e . g. iLock wood et al., 2010; _[Schriiver et al 



20 1 1[ IJones et all \2Q 1 2t iRozanov et all |20 1 2h~ 
4.3 Discussion of CCM results 

We described the impact of NRLSSI and SORCE data, which 
represent the lower and upper boundaries of SSI solar cy- 
cle estimated variations, on the atmosphere and climate as 
depicted in CCM simulations. The NRLSSI reconstructions 
provide the standard data base for simu l ations of the recent 
past a nd future (e.g. ISPARC-CCMVall l2010t iTavlor et al. 
20121) . The atmospheric response with respect to this stan 



dard data set is compared to that derived from a different SSI 
estimate to understand not only the single model responses, 
but also to point out the importance and robustness of so-1905 
lar cycle signals for climate simulations. In particular, it is 
worth mentioning that an enhanced spectral resolution in the 
radiation codes leads to enhanced sensitivity in the response. 
We also described the important role of solar induced ozone 
changes for the amplification of solar effects on atmosphericigio 
composition, circulation and climate. 

Model simulations using the SORCE measurements as 
compared to the NRLSSI model show larger (by a factor 
of two) SW heating and temperature signals. The lower irra- 
diance in the visible range during higher solar activity thani9i5 
at minimum activity in the SORCE measurements does not 
affect the increase in total radiative heating. Recent atmo- 
spheric model simulations with enhanced spectral resolution, 
however, point to the importance of the Chappuis bands in 
lower stratospheric heating. The solar ozone signal derived92o 
when the NRLSSI data are used is positive throughout the 



stratosphere and mesosphere, whereas the sign is reversed in 
the upper stratosphere when the SORCE data are used. Ob- 
servations indicate a reversal of sign in the solar cycle ozone 
response in the upper atmosphere, however, the transition al- 
titude for the sign change is higher than models suggest. This 
result needs to be confirmed by other satellite ozone mea- 
surements as observational evidence of the solar ozone signal 
is still limited. Finally, the tropospheric and surface solar cy- 
cle response has been presented in an ensemble mean and it 
has been highlighted that while these solar induced changes 
are of minor importance for globally averaged temperatures 
there are larger regional responses. 



5 Conclusions 

This paper presents an overview of our present knowledge 
of the impact of solar radiative forcing on the Earth's atmo- 
sphere. It covers the observations and the modelling of the 
solar radiative input as well as the modelling of the Earth's 
atmospheric response. The focus is on satellite-era (i.e. post- 
1970) data for which direct solar irradiance observations are 
routinely available. Special attention is given to the role of 
the UV spectral region, whose small contribution to TSI is 
compensated by a high relative variability with a potentially 
amplified influence on climate through radiative heating and 
ozone photochemistry. There is today clear evidence for 
the impact of solar variability on climate but both its mag- 
nitude and its confidence level are still subject to consider- 
able debate. One major challenge lies in the extraction of 
the weak solar signal from the highly variable atmospheric 
state. Recent progress has been made along two directions. 
The first one is the assessment of the magnitude of secu- 
lar trends in solar radiative forcing through the reconstruc- 
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tion of solar activity on centennial and millenial time scales 
from indirect indices such as cosmogenic isotopes , (e.g . 



Solanki et al.Ll2004tlBard and Frank! 120061: iBeer et al.L 12006; 
Usoskinll2008tlSchmidt et al.Ll2012l) . The second one, which 
we focus on, deals with shorter time scales, and is about re- 
cent solar variability and its impact on the lower and middle 1 
atmosphere. The global physical mechanisms that cause the 
solar irradiance to change in time and ev entually impact cli- 



mate have been well documented, (e.g. lHaigh et al.L I2005 



Haighl I2007L iGrav et al.L [2010t iLean and Woodsl l2010i and 



references therein). More than three decades of SSI observa- 
tions are now available, but they are highly fragmented and 
agree poorly because of the difficulties in making radiometri- 
cally calibrated and stable measurements from space. Merg- 
ing these different observations into one single and homo- 
geneous record is a major and ongoing effort. As a conse- 1990 
quence, several models have been developed for reproduc- 
ing the SSI and its variability. The most successful ones 
are semi-empirical models such as SATIRE, SRPM, and 
COSI, which describe the SSI in terms of contributions com- 
ing from different solar surface magnetic features such as 1995 
sunspots and faculae. Most models nowadays reproduce SSI 
measurements on short term time scales fairly well. How- 
ever, uncertainties in SSI changes still remain on long term 
time scales and in the 220^400 nm band, which is of par- 
ticular interest because of its impact on stratospheric ozone_, 000 
These modelled or observed variations in the SSI are today 
used as inputs to CCM simulations that are capable of prop- 
erly reproducing most aspects of stratospheric heating and 
point to the existence of a significant impact of solar vari- 
ability on climate. However, major uncertainties remain in 
their detailed description, in which nonlinear couplings and 2005 
regional effects can play an important role. 

The main topics discussed in this paper and main conclu- 
sions of our study are: 

- The spectral and temporal coverage of the SSI measure-^ 
ments gradually improved over time, and with the full 
operation of the SORCE mission in April 2004, daily 
observations of the full UV, visible and NIR spectrum 
became available. Unfortunately, this situation is likely 
to end in 2013, when SORCE is anticipated to suc- 
cumb to battery failure. The lack of SSI observationssms 
has led to intensive application of semi-empirical mod- 
els. There remains a considerable issue in assimilating 
SSI observations in such models and in reconstructing 
the SSI prior to the space age. Reconstructions going 
back to the early 20th century can be derived from semi- 
empirical models ba sed on historical gro und-based so- 2020 
lar observations (e.g. Ermolli et al.L 120091) . 



- Recent SSI measurements by SORCE suggest a larger 
(factor 2-6) variability in the UV (200-400 nm) during 
solar cycle 23, which is hard to reconcile with earlier2025 
SSI and TSI observations (e.g. by UARS/SOLSTICE 
and UARS/SUSIM) and with SSI models. However, 



new estimates of possible calibration corrections sug- 
gest that the UV irradiance variation derived from 
SORCE measurements migh t be almost a f actor of two 
weaker than reported earlier (IWoodsl 120121 . 



TSI alone does not adequately describe the solar forc- 
ing on the atmosphere and therefore SSI variations 
have to be taken into account in climate models. For 
many years, the canonical value of the average TSI 
was 1365.4 ± 1.3 Wm~ 2 whereas now the most accu- 
rate, and general l y accepted, value is 1361 ± 0.5 Wm -2 



( Kopp and Leari 201 1 



Sch mutz et all 1201 2) 



There has been steady progress in modelling SSI varia- 
tions, and a number of models are now available that 
can be used as input for climate studies. However, 
the main uncertainty in the models concerns the wave- 
length range 220-400 nm, where the magnitude of the 
variations differs by as much as a factor of three be- 
tween models. This range is of particular interest for cli- 
mate studies. Of the five SSI models discussed, specif- 
ically NRLSSI, SATIRE-S, COSI, SRPM and OAR, 
only one (SRPM) shows a behaviour of the UV and 
visible irradiance qualitatively resembling that of the re- 
cent SORCE/SIM measurements. However, it should be 
noted that the integral of the SSI computed with SRPM 
over the entire spectral range (i.e. the TSI) does not re- 
produce the measured cyclical TSI changes. None of the 
other four models, which are in closer agreement with 
each other, reproduces the peculiar behaviour of the UV 
and visible irradiance observed by SORCE/SIM. 

While there has been major recent progress in better re- 
producing the SSI changes on short term time scales, 
there remains now an important issue in the derivation 
of realistic confidence intervals, both for the observa- 
tions and for the model results. Further common metrics 
should be used for comparing them. 

Within the range of recent SSI values from observa- 
tions and semi-empirical models, the NRLSSI model 
and SORCE observations represent, respectively the 
lower and upper limits in the magnitude of the solar cy- 
cle variation in the UV. 

Results obtained with CCMs show that the observed in- 
crease in UV irradiance at solar activity maximum com- 
pared to minimum leads to an increase in atmospheric 
heating rates and correspondingly an increase in strato- 
spheric temperatures. This direct atmospheric response 
is larger for larger UV forcing, i.e. for SORCE as com- 
pared to NRLSSI, and is sensitive to resolution in the 
radiation code. Currently there is insufficient observa- 
tional evidence to support recent SSI measurements by 
SORCE on the basis of comparisons between climate 
model simulations and atmospheric ozone or temper- 
ature observations in the stratosphere. The larger UV 
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forcing also leads to a larger surface response. The sur- 
face effect is regional and has little influence on globally 
2030 averaged temperatures. 

2085 

- Accurate representation of the spectral nature of the 
incoming radiation, especially at wavelengths below 
320 nm, and therefore ozone photochemical variations 
in the model simulations, is important since these 
2035 changes amplify the atmospheric solar signal. A more 2090 

accurate representation of the Chappuis absorption band 
in a radiation scheme revealed enhanced sensitivity in 
the heating rate response, which again highlighted the 
importance of a well resolved radiation code. 2og5 

2040 A unique aspect of this study is the description of the solar 
terrestrial connection by an interdisciplinary team of solar 
and atmospheric physicists. Progress on this hotly debated 
issue has often been hampered by the fact that limitations 

r J 2100 

on observations or on models are not always properly known 

2045 outside of a given scientific community. For the first time 
a comprehensive comparison and discussion of all relevant 
SSI measurements and models available for climate studies 
is presented, as well as a first investigation of their impacts 
on Earth's climate within a number of different CCMs. These 2105 

2050 results highlight the importance of taking into account in fu- 
ture climate studies SSI variations and their effects on the 
Earths atmosphere. 

Major efforts, however, are still needed in each of the three 
research areas covered in this study. Resolving differences 

2055 in existing long-term time series of solar measurements is a 
major challenge; forthcoming initiatives for merging these 
measurements into single homogenous databases will be ofeno 
invaluable help (see, e.g. http://projects.pmodwrc.ch/solid/). 
Overcoming inaccuracies of current SSI models is another 

2060 important challenge that will benefit from the renewed inter- 
est in the solar radiative output, following the recent solar 
minimum. Finally, the realistic evaluation by climate mod- 2115 
els of the solar impact on the Earths climate is an important 
issue. For example, we focussed on the effect of the solar 

2065 forcing without quantifying the impacts on amplification and 
feedback mechanisms. This should be done in a coordinated, 120 
set of CCM experiments where the treatment of SSI inputs to 
the models are completely specified and results are robustly 
comparable with each other. Then it will be also possible 

2070 to investigate the effects of the top-down feedback and for 
CCMs with an interactive ocean also the bottom-up feedback 2125 
mechanism. 
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